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ABSTRACT:  The  location  of  power  reactors  in  urban  areas 
Increases  the  need  for  containment  of  the  raaximxan  credible 
nuclear  accident.  Such  an  accident  could  possibly  Involve 
a  large  and  rapid  energy  release  in  the  core  region  due  to 
a  nuclear  excursion.  To  contain  the  energy  i*elease  (and 
any  subsequent  release  of  fission  products),  the  Integrity 
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The  work  described  In  this  report  was  performed  under  Tasks 
I  -  IV  of  NOL  Task-285,  NOL  Reactor  Vessel  Containment  Program, 
during  the  period  1956-1964.  The  primary  objective  of  these 
tasks,  to  determine  the  optlmxmi  containment  design  of  nuclear 
reactor  vessels  by  utilizing  the  mechanism  of  gross  plastic 
deformation  as  a  means  of  energy  absorption,  was  accomplished 
along  with  the  successful  achievement  of  a  number  of  secondary 
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CHAPTER  1 
INTRODUCTION 

1.1  Foreword.  The  Reactor  Vessel  Containment  Program  was 
conducted  at  the  U.  S.  Naval  Ordnance  Laboratory,  White  Oak 
and  sponsored  by  the  Atomic  Energy  Commission's  Research  and 
Development  Branch  over  a  period  of  eight  years  (1956-1964). 

The  program  was  a  fundamental  and  comprehensive  research 
effort  directed  to  defining  the  magnitude  of  a  nuclear  excursion 
(in  terms  of  TNT  energy  and  flux)  that  a  primary  nuclear  vessel 
Is  capable  of  containing.  The  results  of  the  work  provide 
relations  that  prescribe  the  design  specifications  of  reactor 
vessels  for  containing  the  energy  releases  due  to  such  excursions. 

1.2  Reactor  Accidents.  In  considering  the  safety  of  any  reactor 
system,  a  large  number  of  potential  accidents  are  examined  In 
the  light  of  their  possible  consequences.  That  accident  with 
the  worst  consequences  for  the  public  safety  and  which  Is  con¬ 
sidered  credible  within  the  bounds  of  a  particular  reactor 
design  Is  often  termed  the  mauclmum  credible  accident  or  MCA 

(as  opposed  to  Incredible,  although  possibly  more  hazardous, 
due  to  the  Inherent  reactor  characteristics  or  mechanisms). 

This  report  considers  a  type  of  accident  that  Is  sometimes 
considered  credible,  namely  a  large  and  rapid  release  of  energy 
at  the  reactor  core  due  to  an  uncontrolled  nuclear  excursion. 
Whether  or  not  such  an  accident  Is  the  MCA  depends  upon  the 
particular  reactor  design. 

The  excursion  energy  generated  at  the  reactor  core  Is 
frequently  expressed  In  terms  of  megawatt-seconds;  however, 
because  of  similarities  In  rate  of  energy  release,  this  acci¬ 
dent  may  also  be  expressed  In  pounds  of  TNT  particularly  where 
destructiveness  and  containment  Integrity  are  basic  questions. 
Bethe  and  Talt  In  reference  (a),  for  example,  were  among  the 
first  to  characterize  nuclear  accidents  In  terms  of  TNT. 
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1.3  Reactor  Containment.  W.  B.  Cottrell  presents  a  somewhat 
detailed  containment  philosophy  In  reference  (b);  only  a  brief 
treatment  will  be  given  here.  Containment  Is  a  provision  to 
control  (limit)  the  public  exposure  to  radiation  following  any 
reactor  accident  to  acceptable  levels.  Early  In  nuclear  develop¬ 
ment  such  control  was  often  aided  In  the  United  States  through 
locating  reactors  In  huge  and  generally  remote  areas  such  as 

Oak  Ridge  and  the  desert  regions  of  Washington  and  Idaho.  In 
these  and  similar  Instances,  the  distance  from  each  reactor  to 
the  site  boundary  was  intended  to  permit  adequate  control  in 
the  event  of  an  accident. 

Several  needs  have  evolved  to  radically  change  this  picture. 
The  most  significant  of  these  Is  that  location  In  wilderness 
areas  is  not  consistent  with  today's  electric  power  requirements 
for  our  industrial  complexes,  since  these  requirements  generally 
Imply  location  near  or  even  In  the  midst  of  heavily  populated 
areas.  When  central  power  station  requirements  are  coupled  with 
the  growing  use  of  atomic  plants  for  marine  propulsion  and  for 
research  purposes,  reactor  containment  Is  brought  directly  to 
the  ’’front  door”  of  the  general  public.  Consequently,  contain¬ 
ment  systems  consisting  of  combinations  of  structural  barriers 
and  mechanical  auxlUarles  must  be  employed  to  restrict  the  pos¬ 
sible  spread  of  radioactive  fission  products  that  might  be 
released  from  the  core  by  an  accident.  One  such  barrier  might 
be  considered  to  be  the  reactor  primary  vessel  If  Its  Integrity 
could  be  guaranteed. 

1.4  Scope  and  Significance.  The  Reactor  Vessel  Containment 
Program  was  directed  to  characterizing  the  excursion  containment 
potential  of  nuclear  facilities  In  general  and  reactor  primary 
pressure  vessels  In  particular;  this  was  the  broad  purpose. 

The  ultimate  objective  was  to  devise  general  laws  for  containment 
of  the  excursion  within  the  reactor  vessel  and  to  express  these 
laws  In  a  fundamental  design  format  that  would  permit  Immediate 
and  fruitful  use  by  containment  engineers. 
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The  theoretical  and  experimental  studies  of  the  program 
were  extensive  In  spectrum  and  depth.  Since  the  excursion 
loading  and  vessel  response  parameters  were  so  numerous  amd 
so  complex,  Salnt-Venant 's  semi-inverse  method  of  classical 
mechanics  was  Invoked  to  generate  an  embryonic  functional  form 
of  the  desired  solution,  ITie  necessary  boundary  conditions 
were  obtained  from  over  100  experiments  employing  explosives 
and  propellants  In  scaled  model  cylinders  simulating  reactor 
vessels.  What  follows  Is  an  unfolding  of  the  rationale  and 
the  broad  spectrum  of  tools  -  both  theoretical  and  experimental 
that  were  employed  to  characterize  the  parameters  of  containment 
and  to  formulate  excursion  or  EXPLOSION  CONTAINMENT  LAWS  FOR 
NUCLEAR  REACTOR  VESSELS. 
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CHAPTER  2 

ANCILLARY  PROGRAMS 

The  principal  mission  of  the  Reactor  Vessel  Containment 
Program,  as  originally  planned,  was  to  conduct  basic  studies 
of  a  general  character  directed  to  containment  per  se.  Prom 
time  to  time,  however,  problems  with  particular  reactors  arose; 
these  problems  and  their  solutions  are  recoxmted  briefly  In 
the  following  paragraphs  because  they  have  provided  Information 
basic  to  the  development  of  Containment  Law  rationale. 

2.1  Preliminary  Stvtdy  of  Enrico  Fermi  Plant.  In  1957,  upon 
request  of  the  Atomic  Energy  Commission,  program  effort  was 
diverted  temporarily  to  a  preliminary  containment  study  of  the 
Enrico  Perml  Atomic  Power  Plant.  The  results  of  this  study 
were  reported  by  E.  M.  Plsher  emd  W.  R.  Wise  In  reference  (c); 
they  found  that : 

(a)  An  excursion  releasing  energy  equivalent  In  magnitude 
and  violence  to  1000  pounds  of  TNT  detonated  at  the  core  would 
not  develop  shook  waves  of  sufficient  strength  or  missiles  with 
sufficient  energy  to  threaten  the  Integrity  of  the  secondary 
shield. 

(b)  The  response  of  the  rotating  shield  plx;g  to  the  shock 
would  be  a  vertical  Jump  not  exceeding  two  feet. 

(c)  But  the  Internal  blast  pressure,  the  quasi  static 
pressure  generated  by  the  gases  released  from  the  explosive 
charge,  would  be  the  predominate  mechanism  that  could  propel 
the  massive  shield  plug  throxigh  the  roof  of  the  containment 
building. 

2.2  Energy  Partition  Study.  The  Perml  plug-Jvmp  hazard  reported 
In  reference  (c)  raised  serious  questions  relating  to  disposition 
of  the  energy  during  a  nuclear  excursion.  In  consequence,  upon 
AEC  request  to  the  Naval  Ordnance  Laboratory  In  1958,  It  was 
decided: 
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To  determine  the  partition  of  meohanloal 
and  non -mechanical  energy  resulting  from 
the  detonation  of  water- cased  and  sodium- 
cased  explosives  surroundtd  by  air  and 
confined  within  a  closed  plxig-fltted  vessel. 

The  Initial  efforts  made  on  this  task,  reported  by  Wise 
In  reference  (d)  and  by  J.  P.  Proctor  In  detail  In  reference 
(e),  were  exploratory  In  nature  and  served  principally  to 
establish  problem  approach  and  the  governing  parameters  of 
the  problem.  Parameters  such  as  charge  weight,  mass-per- 
frontal  area  of  the  shield  plvig,  llquld-to-alr  ratio,  and 
liquid  temperature  were  examined  In  an  exploratory  apparatus 
that  represented  a  crude  model  of  the  Enrico  Perml  Reactor, 

Three  Important  conclusions  derived  from  this  work  were: 

(a)  The  presence  and  temperature  of  a  liquid  casing 
surrounding  the  charge  greatly  affected  energy  partition: 
Increasing  quantities  of  liquid  heat  sink  decreased  the  plug- 
jump  and  Increasing  temperatures  Increased  the  Jump. 

(b)  A  refined  experimental  apparatus  would  be  required  for 
accurate  characterization  of  the  principal  parameters. 

(c)  The  explosion  containment  potential  of  the  Fermi 
composite  envelope  (reactor  vessel,  graphite  shield,  and 
primary  shield  tank)  must  be  established  If  a  meaningful  plug- 
response  solution  were  to  be  achieved, 

2.3  Extensive  Perml  Shield  Plug  Study.  The  previously  dis¬ 
cussed  energy  partition  studies,  comprising  destructive 
mechanisms  basic  to  reactor  excursions  generally,  were  con¬ 
tinued  until  October,  1959.  At  that  time.  Atomic  Power 
Development  Associates,  Inc.  (APDA)  -  the  designers  of  the 
Perml  Plant  -  requested  the  AEG  to  assist  In  resolving  the 
plug- Jump  question  as  quickly  as  possible.  Specifically, 

APDA  reconunended  that  MOL  conduct  a  program  addressed  to  the 
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excursion  motion  of  the  Perml  shield  plug  per  se.  Accordingly, 
upon  AEG  request  to  NOL,  the  general  energy  partition  studies 
were  discontinued,  and  the  entire  manpower  and  funds  of  the  NOL 
Containment  Program  were  directed  to  the  following  particularized 
objective : 

To  determine  the  response  of  the  Perml 
shield  plug  to  a  postulated  nuclear 
excursion  equivalent  In  energy  release 
and  rate  of  release,  l.e,,  equivalent  In 
violence  and  destructiveness,  to  the 
detonation  at  standard  conditions 
(Temp:  77®P>  Press:  atmospheric)  of  a 
center-initiated,  compact  cylindrical 
charge  of  1000  pounds  of  TNT  located  at 
the  reactor  core. 

Preliminary  results  of  this  extensive  containment  study 
on  the  Perml  Reactor  (1959-1962)  were  reported  by  Wise  In 
reference  (f),  by  Wise,  Proctor,  and  Walker  In  reference  (g), 
aind  by  H.  B.  Beneflel  In  reference  (h).  The  final  and  compre¬ 
hensive  results  of  the  entire  Pennl  program  were  reported  by 
Wise,  Proctor,  and  Walker  In  reference  (l);  some  of  the  major 
conclusions  were 

(a)  The  maximum  Jump  of  the  l43-ton  shield  plx:g  resulting 
from  the  1000- lb  TNT  accident  would  be  of  the  order  of,  but  not 
exceed,  102  feet.  Thvis  the  plug  could  possibly  be  propelled 
throxigh  the  roof  of  the  containment  bxilldlng. 

(b)  Containment  of  the  postulated  accident  within  the 
reactor  vessel  would  be  marginal  for  no  loss  of  sodium  coolant. 
Por  major  losses  of  coolant,  however,  containment  would  be 
conservative.  In  any  event,  the  gross  containment  Integrity 
of  the  secondary-shield  wall  would  not  be  violated. 
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(c)  Release  of  the  postulated  energy  at  a  rate  slower 
than  the  detonation  rate  of  WT  would  yield  a  plug  Jump  less 
than  that  occuirlng  from  the  TUT  accident. 

(d)  The  graphite  shield,  due  to  Its  Inertia,  Increases 
the  TIfT  explosion  containment  potential  of  the  composite 
(reactor  vessel -primary  shield  tank)  envelope  by  a  factor  of 
two  or  greater, 

(e)  The  sodl\im,  blanket,  thermal  shield,  and  other 
reactor  vessel  Internals  possess  a  significantly  large  capacity 
for  absorbing  heat  from  the  TNT  explosion  gases  during  the 
short  time  that  would  be  required  for  the  shield  pl\]g  to  escape 
the  composite  envelope. 

The  Enrico  Fermi  Reactor  has  been  modified  and  fitted 
with  an  energy-absorption  mechanism  In  the  machinery  done  (see 
reference  (j))  to  restrict  any  upward  plug  motion  resulting 
from  an  accident.  If  the  absorber  possesses  the  potential 
specified  In  reference  (j),  motion  of  the  shield  plug  as  a 
result  of  the  postulated  accident  studied  by  NOL  will  not 
constitute  a  serious  hazard. 
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CHAPTER  3 

THE  PRINCIPAL  PROGRAM 

The  principal  containment  program  consisted  of  several 
related  studies  conducted  over  a  span  of  about  eight  years. 

The  basic  rationale  of  these  studies  and  their  synthesis  follow, 

3.1  Long-Term  Objective.  The  long-term  technical  objective 
of  the  MOL  Reactor  Vessel  Containment  Program  was  to  achieve 
the  optimum  containment  design  of  nuclear  reactor  vessels 
through  gross  plastic  deformation  as  a  means  of  absorbing 
excursion  energy.  The  steps  of  the  objective  were: 

(a)  To  simulate  nuclear  excursions  occurring  In  full- 
size  reactor  vessels  through  the  use  of  propellants  and  explo¬ 
sives  In  model  vessels 

(b)  To  Investigate  reactor  vessel  response  In  the  gross 
plastic  domain  with  material,  size,  shape,  and  constraint  as 
parameters 

(c)  To  formulate  general  laws  for  containment  of  explosion- 
type  energy  releases  within  the  reactor  vessel  and  to  express 
these  laws  In  a  fundamental  design  format  that  will  permit 
Immediate  use  by  containment  engineers. 

Progress  on  the  long-term  objective  was  first  reported 
formally  (1958)  by  Wise  In  reference  (k):  he  derived  the 
General  Strain  Energy  Equation  of  Equilibrium  for  a  Vessel 
Subjected  to  Internal  Dynamic  Loading,  conducted  basic  experi¬ 
ments  with  solid  propellant  In  model  reactor  vessels,  and 
concluded  that  strain  energy  absorption  was  a  significant 
containment  parameter, 

3.2  General  Strain  Energy  Equation.  Consider  a  closed  right- 
cylindrical  vessel,  e.g.,  a  reactor  pressure  vessel,  subjected 
to  an  accidental  excursion  that  causes  the  vessel  wall  to  flow 
rapidly  and  gross -plastically.  If,  during  the  dilation. 
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(a)  the  pressure  vector  at  every  differential  surface 
of  the  Internal  wall,  and 

(b)  the  strain  vector  of  every  differential  surface  of 
the  external  wall 

are  known,  what  Is  the  strain  energy  correlative  with  containment? 

Work  done  upon  the  Internal  wall  of  the  vessel  must  be 
Instantaneously  conserved  In  the  Immediate  mechanical  forms  of 
strain,  kinetic,  compression,  and  potential  (position)  energy, 

A  consequence  of  this  conservation  Is  the  General  Strain  Energy 
Equation  of  Dynamic  Equilibrium  that  was  derived  In  reference  (k) . 
Plastic  straining  of  the  vessel  wall  Is  accompanied  by  the 
generation  of  heat  within  the  wall  which  may  or  may  not  be 
dissipated  from  the  wall  depending  upon  the  temperature  of  Its 
environment.  The  heat  of  elastic  and  plastic  strain  Is  accounted 
for  Inherently  In  the  strain-energy  function.  Although  the 
magnitudes  of  the  surface,  strain.  Inertia,  compression,  and 
position  energies  are  functions  of  the  temperature  of  the  wall, 
the  validity  of  the  equation  Is  Independent  of  temperature  and 
heat  flux  within  the  range  for  which  utilization  of  the  wall  as 
a  pressure  vessel  remains  practicable. 

The  General  Strain  Energy  Equation  Is  greatly  simplified 
for  the  case  of  the  homogeneous  and  Isotropic,  right-circular 
cylinder  closed  with  rigid,  radial  constraints  at  the  ends. 

The  conditions  of  circularity,  cyllndrlcallty,  and  rightness 
of  the  vessel  provide  that  the  spatial  limits  of  Integration 
be  constant.  If  the  loading  Is  rotatlonally  symmetric,  the 
structural  response  of  the  homogeneous.  Isotropic  wall  will  be 
rotatlonally  synmetrlc  -  a  finding  well  established  by  experlMuats 
reported  herein  and  In  reference  (k). 

Since  the  rationale  of  the  General  Strain  Energy  Equation 
Is  so  germane  to  this  report  and  to  containment  generally.  Its 
complete  derivation  and  correlative  material  are  reproduced  In 
Appendix  A. 
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3.3  Strain  Energy  and  Containment.  The  mathematical  and 
physical  formulation  of  the  Strain  Energy  Equation  Is  not  a 
theory  designed  for  experimental  confirmation;  but  rather  a 
principal  analytic  tool  through  which  experimental  work  can 
be  used  to  assess  the  effectiveness  of  strain-energy  absorption 
as  a  means  for  Increasing  the  containment  potential  of  reactor 
vessels.  This  assessment  can  be  expressed  as  the  partition  of 
mechanical  (strain  energy)  and  non-mechanical  energy  resulting 
from  excursion- type  energy  releases  In  simulated  reactor  vessels. 
With  the  total  energy  Input  known,  the  General  Strain  Energy 
Equation  provides  the  means  for  obtaining  the  desired  partition 
without  further  knowledge  as  to  the  disposition  of  non-mechanical 
energy. 

The  model  vessels  subjected  to  simulated  excursion  loading 
In  reference  (k)  were  right-circular  cylinders  with  length  equal 
to  two  diameters,  and  the  material  was  type  304  stainless  steel. 
The  propellant  employed  In  these  basic  experiments  produced 
loadings  generated  In  milliseconds.  The  results  Indicated  that 
the  energy  containable  In  the  plastic  range  exceeded  that  In 
the  elastic  range  by  an  "elastic-plastic"  factor  of  about  eight. 
Prom  the  strain  energy  equation,  the  deformation  energy  was 
found  to  be  about  one-tenth  of  the  total  energy  released. 

For  microseconds  loading  mechanisms  such  as  the  TNT  shook 
wave.  It  was  found  later  that  the  elastic-plastic  factor  exceeded 
that  for  the  propellant  and  that  the  deformation  energy  was  of 
the  order  of  one-half  the  total  energy  released*.  Thus,  the 
valuable  containment  role  of  plastic  strain  energy  became 
manifestly  clear,  and  It  was  concluded  that 

*  The  energy  released  by  a  TNT  explosion  varies  from  1050  to 
3620  calorles/gm  of  TNT  depending  upon  the  availability  of 
oxygen  from  the  surrounding  atmosphere.  For  purposes  of  this 
report.  It  Is  assvmed  thct  the  amount  of  gas  confined  within 
the  reactor  primary  vessel  Is  small  and  Is  generally  Inert  or 
low  In  oxygen  content.  Therefore,  the  energy  released  by  a 
TNT  explosion  within  a  reactor  vessel  would  be  equivalent 
essentially  to  the  heat  of  detonation  -  1050  cal/gm. 


10 


NOLTR  63-140 


Strain  Energy  Absorption,  achieved  through 
the  medium  of  large  plastic  deformation. 

Is  a  means  for  substantially  Improving  the 
containment  potential  of  nuclear  reactor 
vessels,  ancillary  container  paraphernalia, 
and  building  enclosures. 

3.4  Need  for  Containment  Laws.  Philosophically,  containment 
must  be  provided  to  reduce  the  consequence  of  a  reactor  accident 
so  that  the  risk  of  reactor  operation  to  the  public  can  be 
decreased  to  acceptable  levels.  Since  power  reactors  located 
In  urban  locales  cannot  provide  safety  throtigh  Isolation,  they 
must  rely  upon  ’’mechanical”  containment  devices.  In  the  Enrico 
Perml  Plant,  for  example,  there  are  four  such  barriers:  the 
reactor  vessel,  the  primary  shield  tank,  the  massive  secondary 
shield,  and  the  so-called  containment  building.  Of  these, 
however,  the  reactor  vessel  possesses  the  greatest  strain 
energy  absorption  potential  by  far  -  hence.  It  constitutes  a 
significant  containment  barrier  as  shown  In  reference  (l). 

The  reactor  vessel.  In  general,  has  great  mass  and  strength 
dictated  by  operating  design  requirements,  and  these  features 
constitute  desirable  containment  attributes.  It  should  be 
clearly  established  at  the  outset,  however,  that  the  optimum 
operating  design  of  nuclear  reactors  does  not  In  general  con¬ 
stitute  optlm\mi  containment  design.  The  operating  design 
characteristics  of  a  reactor  Inherently  restrict  the  structure 
to  stresses  within  the  elastic  limit  of  the  material.  Accident 
containment  potential,  however.  Is  greatly  enhanced  through 
gross  plastic  deformation.  Although  such  deformation  would 
render  the  reactor  Inoperable,  a  vastly  more  Important  objec¬ 
tive  will  have  been  achieved  If  the  excursion  Is  contained  or 
attenuated  and  a  major  catastrophe  averted.  This  Is  true 
Irrespective  of  the  overall  Integrity  of  the  primary  system, 
since  energy  absorbed  via  gross  plastic  deformation  Is  \mavallable 
for  further  destruction. 
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The  reactor  vessel  can  and  should  be  a  vital.  If  not 
paramoxmt,  containment  device  for  all  power  reactors.  Even 
without  planned  containment  design,  some  reactor  vessels 
possess  a  large  potential  for  absorbing  excursion  energy; 
with  planned  design,  the  containment  potential  of  reactors 
now  on  the  drawing  board  could  perhaps  be  greatly  Increased. 

In  many  cases.  If  excursions  could  be  completely  contained 
within  the  reactor  vessel,  the  requirements  for  additional 
barriers  could  be  reduced  and  the  total  outlay  for  safety 
sharply  curtailed. 

The  preceding  discussion  on  strain  energy  absorption  and 
reactor  containment  potential  is  largely  qualitative  or  seml- 
quantitatlve.  But  engineers  require  highly  quantitative 
Information  for  reactor  vessel  containment  design;  they  must 
know  the  dimensions,  shape,  and  constraints  correlative  with 
containing  postulated  excxirsions.  In  other  words,  they  require 
EXPLOSION  CONTAINMENT  LAWS  FOR  REACTOR  VESSELS  -  laws  expressed 
In  terms  of  conventional  materials  and  properties , 
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CHAPTER  4 

SIMULATION  OP  NUCLEAR  EXCURSIONS 

The  magnitude  and  flux  of  excursion  energy  and  attendant 
gases  determine  the  loading  at  the  Internal  wall  of  the  vessel. 
The  stress-strain  state  In  the  wall  and  the  ultimate  response 
of  the  structure  are  complex  functions  of  the  loading,  functions 
obtainable  only  through  experimental  means  and  feasibly  so  only 
through  experiments  In  scaled  models.  Clearly,  the  character 
of  excursion  loading  and  the  means  for  simulating  such  loading 
had  to  be  considered  further  as  a  first  step  In  the  formulation 
of  Containment  Laws.  This  was  accomplished  through  the  use  of 
conventional  explosives  and  through  further  study  of  previous 
propellant  experiments. 

4.1  TNT  Rationale.  The  rationale  for  employing  TNT  (or  a 
similar  conventional  explosive  like  pentollte)  to  simulate 
nuclear  accidents  can  be  readily  expressed:  Nuclear  excursions 
cannot  be  micro-scaled;  TNT  explosions  can,  and, furthermore, 
many  of  their  scaling  properties  are  well  known.  This  makes 
available  a  very  Important  tool  that  Is  often  the  sine  qua  non 
for  complex  multi -parameter  problems :  experiments  In  small 
models.  Although  the  character  of  the  energy  released  during 
a  postulated  nuclear  excursion  may  not  be  Identical  to  that 
of  the  TNT  explosion,  the  violence  factors  of  the  two  mechanisms 
are  often  sufficiently  alike  to  permit  use  of  the  TNT  model  as 
a  reasonable  upper-bound  simulant  (reference  (a)). 

For  reactors,  postulated  credible  accidents  cover  a  wide 
spectrvnn  In  both  magnitude  and  flux  of  the  energy  released; 
these  range  from  small  releases  to  thousands  of  megawatt-seconds 
occurring  In  times  from  a  few  microseconds  to  seconds  or  more. 

In  this  respect.  It  Is  singularly  Important  to  note  that  conven¬ 
tional  high  explosives,  like  TNT  and  pentollte,  also  produce 
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flux-rate  mechanlBms  that  cover  a  broad  apectrum:  the  ahock 
wave  of  the  conventional  exploalve  la  a  ralcroaeconda  mechanlam, 
and  the  quaal-Btatlc  ’’equilibrium”  preaaure  generated  by  the 
releaaed  gaaea  occurs  In  the  domain  of  mllllaeconds  and  seconds. 
Since  the  fundamental  characteristics  of  both  mechanisms  are 
known,  their  application  to  btuslc  studies  Is  feasible.  This 
was  demonstrated  In  the  extensive  structural  response  Investiga¬ 
tions  of  the  Enrico  Fermi  Reactor  previously  referred  to  In 
Chapter  2. 

At  the  energy  density  of  1050  calories  per  gram  -  a 
nominal  value  for  the  heat  of  detonation  of  TNT  -  a  pound  of 
TNT  would  release  two  megawatt-seconds  of  energy.  Since  the 
TTfT  energy  flux  constitutes  a  reasonable  upperbound  for  nuclear 
excursions,  there  exists  the  practice  of  expressing  large 
excursion  accidents  In  pounds  of  TNT  for  purposes  of  safety 
assessment.  That  this  practice  Is  profitable  In  containment 
studies  Is  one  of  the  principal  theses  of  this  report. 

4.2  The  TNT  Explosion.  Consistent  with  the  theme  of  Interest, 
It  Is  assmned  that  the  simulated  excursion  would  be  contained 
within  the  reactor  vessel.  It  Is  further  assumed  that  normal 
operating  pressures  and  temperatures  exist  everywhere  In  the 
reactor  vessel  and  Its  environs  Immediately  prior  to  detonation. 
The  TNT  charge  Is  located  at  the  core;  It  has  a  weight  density 
of  about  1.5  gm/cm3,  it  occupies  a  volume  of  the  order  of  the 
core  volume,  and  Its  temperature  and  pressure  are  25°C  emd 
atmospheric. 

Upon  detonation,  the  detonation  front  would  propagate 
through  the  charge  with  a  velocity  of  about  6OOO  metei»B/Bec 
and  pressures  behind  the  front  estimated  to  be  of  the  order  of 
200  kllobars.  In  consequence,  the  solid  explosive  would  be 
transfozmied  almost  dlso(xitlnuously  In  time  to  a  gaseous  fireball 
at  virtually  the  same  volume  with  temperatures  and  pressures  of 
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the  order  of  3000°C  and  50,000  atmospheres  (see  Cole, 
reference  (l)).  The  energy  density  and  gas  produced  per 
gram  of  TNT  would  be  about  I050  calories  and  0.033  gram 
mole.  As  a  consequence  of  the  immense  energy  flux,  two 
destructive  mechanisms  would  be  produced:  the  shock  wave 
and  the  internal  blast  pressure,  i.e.,  the  quasi'Static 
(or  equilibrivm)  pressure  generated  by  gases  released  from 
the  explosive  charge. 

4.3  Shock  Wave.  The  shock  wave,  which  initially  contains 
approximately  one-half  of  the  energy  released,  would  propa¬ 
gate  through  the  media  surrounding  the  core  and  strike  the 
internal  wall  of  the  reactor  vessel.  In  response  to  the 
delivered  impulse,  the  wall  would  be  given  a  virtually 
instantaneous  velocity.  Hils  is  true,  since  the  response 
time  of  the  wall  would  be  large  relative  to  the  significant 
duration  of  the  shock  wave.  Portions  of  the  wall  nearest  the 
charge  would  be  given  the  greatest  velocity  and  undergo  the 
greatest  deformation.  The  inertia  and  tensile  strength  of 
the  reactor  vessel  would  resist  the  impulsive  force  during 
acceleration  of  the  wall.  Subsequently,  the  tensile  constraints 
would  bring  the  wall  to  rest. 

In  general,  shock  wave  attenuation  in  elastic  media  is 
less  than  in  inelastic,  plastic,  or  relatively  compressible 
media.  Attenuation  of  the  shock  would  be  greater  in  the 
reactor  liquid  coolant  than  in  steel,  and  in  air  or  similar 
gases  the  attenuation  would  be  considerably  greater  than  in 
the  liquid  or  steel.  Thus,  it  is  seen  that  if  the  coolemt 
were  lost  from  the  reactor  vessel,  the  shock  attentuation 
potential  of  the  system  would  be  increased,  emd  the  shock 
strength  at  the  vessel  wall  would  be  decreased  with  a  conse¬ 
quent  decrease  in  deformation.  Clearly,  then,  the  deformation 
volvime  of  the  reactor  vessel  would  be  a  function  of  reactor 
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coolant  content:  the  largest  defonatlon  would  be  associated 
with  the  case  of  no  coolant  loss  and  conrersely  for  the  oo«- 
plete  loss  of  coolant. 

4.4  Internal  Blast  Pressure.  Following  the  shook  ware  would 
be  the  explosion  gases «  which  Initially  contain  the  remaining 
one-half  of  the  released  energy «  and  these  gases  would  be  con¬ 
fined  In  the  reactor  vessel.  The  resultant  Internal  blast 
pressure,  a  (|uasl-statlo  pressure  as  opposed  to  the  short- 
duration  pressure  of  the  shock  wave,  would  be  a  significant 
function  of  two  variables:  the  volume  Into  which  the  explosion 
gases  flow  and  the  capacity  (quality  and  quantity)  of  the  heat 
sinks  available  to  these  gases.  Distortion  of  the  vessel  wall 
from  the  shock  wave  would  precede  arrival  of  the  blast  pressure 
at  the  wall,  so  that  the  expansion  volume  available  to  the 
explosion  gases  would  be  the  Siim  of  the  Initial  gas  volixme, 

the  void  volume  created  by  coolant  losses,  and  the  deformation 
volume  produced  by  the  shock. 

If  a  major  portion  of  llqxild  coolant  were  retained  In  the 
vessel,  the  coolant  would  constitute  a  significant  heat  sink. 

If  a  major  portion  or  all  of  the  coolant  were  lost  from  the 
vessel,  the  complex  metal  structures  that  generally  surround 
the  core  would  be  Important  heat  sinks.  In  any  event,  these 
heat  sinks  would  considerably  reduce  the  Internal  blast  pressure. 
Through  expansion  and  heat  losses,  the  blast  pressure  would  be 
generally  an  order  of  magnitude  or  more  lower  than  the  shock 
pressure,  and  would  not  be  sufficient  to  further  distort  the 
vessel  wall.  This  Is  shown  In  reference  (l)  and  discussed 
further  In  section  7.6. 

4.5  Propellant  Loading.  Again  It  Is  assumed  that  ccmtalnment 
occurs  and  that  normal  operating  pressuires  and  temperatures 
exist  everywhere  In  the  reactor  and  Its  environs  Immediately 
prior  to  Initiation.  The  solid  propellant  Is  at  the  core;  Its 
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weight  and  energy  densities,  temperature  and  pressure  are  like 
those  for  the  TNT  case  Just  discussed.  Upon  Initiation  and 
burning,  however,  the  character  of  the  loading  produced  wovild 
be  substantially  different  from  that  of  TNT,  since  the  energy 
release  time  would  be  of  the  order  of  milliseconds  and  seconds. 

No  shock  wave  would  be  produced;  consequently,  all  of  the 
released  energy  would  appear  initially  in  the  product  gases. 

As  the  propellant  continued  to  burn,  a  stable  and  uniform  build¬ 
up  of  pressure  would  occur  (see  reference  (k)  and  section  9.7). 
Like  the  internal  blast  pressure  of  TNT,  the  propellant  loading 
would  be  a  significant  function  of  the  expansion  volume  and  the 
capacity  of  the  available  heat  sinks  -  both  of  which  would  act 
to  decrease  pressure.  Since  the  propellant  does  not  generate 
a  shock  wave  that  is  immediately  available  for  vessel  deforma¬ 
tion,  energy  and  gas  release  must  be  large  enough  to  overcome 
the  negating  effects  of  the  heat  sinks  and  expansion  volvime  before 
a  sufficient  pressure  can  be  developed  to  rupture  the  vessel. 

This  pressure  would  correspond  nearly  to  the  static  rupture 
strength  of  the  vessel.  It  is  noted,  therefore,  that,  for 
marginal  containment,  the  post-accident  pressure  in  the  reactor 
vessel  resulting  from  a  propellant  loading  would  exceed  con¬ 
siderably  that  from  a  TNT  explosion  (internal  blast  pressure). 

For  marginal  containment,  the  deformation  pattern  would  be 
similar  to  that  for  TNT,  and  hence,  the  expansion  volvimes  would 
be  of  the  same  order. 

4.6  TNT-Propellant  Efficiency.  The  heat  sinks  in  the  reactor 
vessel  would  attenuate  the  propellant  loading  to  a  much  greater 
degree  than  for  the  TNT  case.  This  follows  since  the  exposure 
time  would  be  large  for  all  of  the  propellant  energy  released, 
whereas  for  TNT,  50  per  cent  of  the  energy  released  is  in  the 
shock  wave  delivered  to  the  vessel  wall,  and  this  is  independent 
of  the  capacity  of  the  heat  sinks.  Indeed,  this  time -dependent 
absorption  of  energy  is  an  overriding  containment  phenomenon 
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from  which  the  following  general  conclusion  can  be  drawn: 
the  vessel  deformation  efficiency  of  TNT  is  greater  than  that 
of  slower  energy  releases;  conversely,  the  reactor  containment 
potential  for  energy  flxuces  less  than  that  of  TOT  is  greater 
than  for  the  TMT  flux. 

The  Influence  upon  containment  of  D’Alembert  (or  Inertial) 
constraints  is  another  point  of  important  difference  between 
propellant  and  TNT  shock  loading.  Very  massive  components 
like  the  143-ton  Fermi  shield  plug  would  be  negligibly  responsive 
to  the  shock,  but  significantly  so  to  propellant  loading.  The 
same  rationale  applies  to  the  massive  graphite  blocks  surround¬ 
ing  the  Fermi  reactor  and  the  large  body  of  water  surroxmdlng 
Oak  Ridge’s  High  Flux  Isotope  Reactor:  the  graphite  and  water 
would  enhance  the  reactor  containment  potential  for  shock  loading 
but  make  almost  no  contribution  to  containment  for  the  slower 
propellant  loading. 

In  summary,  it  is  seen  that  energy  fl\ix  constitutes  an 
extremely  significant  excursion  containment  parameter.  Energy 
release  mechanisms  of  microseconds  and  milliseconds  duration 
have  markedly  different  complexions  and  present  greatly  different 
containment  problems.  And  with  respect  to  "efficiency": 

(a)  the  deformation  efficiency  of  TOT  (the  capacity  for 
doing  work  on  constraints  that  are  largely  non-lnertial)  is 
greater  than  that  for  slower  energy  releases 

(b)  the  D’Alembert  efficiency  of  TNT  (the  capacity  for 
doing  work  on  constraints  that  ai?e  largely  inertial)  is  less 
than  that  for  slower  energy  releases. 

4.7  50/50  Pentolite.  Mixtures  of  PETO  (C5H8O12N4)  and  TNT 

(C7H5O6M3)  are  called  pentolites,  the  most  common  being  a 
50/50  mixture.  All  of  the  previous  discussions  on  the  behavior 
of  TNT  applies  to  pentolite,  and  pentolite  possesses  additional 
advantages  as  well.  In  view  of  the  ease  with  which  it  is 
detonated  and  its  established  replicability  in  small-scale  tests. 
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50/50  pentolite  was  chosen  as  the  most  satisfactory  explosive 
charge  to  employ  in  the  containment  pregram.  The  pentolite/nTP 
ratios  of  energy  density  and  release  of  gas  products  per  unit 
of  reactant  are  1.02  and  1.17^  respectively;  hence  the  use  of 
pentolite  in  containment  studies  is  consistent  with  reasonable 
upper-bound  criteria. 
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CHAPTER  5 

METHOD  OP  SOLUTION 


The  solution  to  the  multi-faceted  containment  problem 
required,  first  of  all,  that  a  Judicious  choice  of  the 
principal  parameters  be  identified  and,  at  least,  partially 
characterized . 

5.1  Principal  Containment  Parameters.  In  review  of  the 
loading  mechanisms  discussed  in  Chapter  4,  it  is  seen  for 
the  reactor  vessel  per  se  that  the  principal  containment 
parameters  can  be  taken  as : 

(a)  the  reactor  vessel  material,  size,  shape,  and  con¬ 
straints 

(b)  the  magnitude  and  duration  of  the  shock  wave  striking 
the  vessel  wall 

(c)  the  inertia  and  dynamic  tensile  strength  of  the 
vessel  wall 

(d)  the  heat-sink  capacity  of  the  internals  of  the 
reactor  vessel,  including  the  liquid  coolant. 

These  basic  quantities  are  Innate  to  excursion  containment. 

They  are  not  all  true  parameters,  however;  some  are  pseudo 
parameters,  consisting  of  a  fixed  element  and  an  element  that 
variee ; 

(a)  The  shock  wave  strength  is  a  function  of  the  density 
of  the  media  surrounding  the  core;  the  density  of  the  steel 
internals  stays  fixed,  but  that  of  the  remainder  varies  with 
coolant  content. 

(b)  The  inertia  of  the  vessel  wall  is  a  function  of  the 
wall  mass  and  acceleration;  the  mass  remains  fixed,  but  the 
acceleration  varies  with  the  strength  of  the  shock  wave  and 
the  dynamic  tensile  strength  of  the  reactor  vessel. 
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(c)  The  dynamic  tensile  strength  of  the  vessel  nail 
depends  upon  the  static  strength  of  the  wall,  gross  strain, 
and  the  rate  of  strain;  the  static  strength  Is  fixed,  but 
the  strain  and  rate  of  strain  vary  with  the  strength  of  the 
shock  wave. 

(d)  The  total  heat-sink  capacity  of  the  reactor  vessel 
Internals  consists  of  the  combined  sinks  afforded  by  the 
reactor  coolant  and  the  metal  structures  that  generally 
surround  the  core;  the  capacity  remains  constant  with  respect 
to  the  metal,  but  varies  with  coolant  content. 

5.2  Fundamental  Questions.  Attempts  to  quantify  and  ssmthe- 
slze  the  above  parameters  raised  some  Interesting  questions. 

(a)  What  are  the  shock  attenuation  properties  of  the 
metal-coolant-gas  medium  that  surrounds  the  core?;  how  are 
these  properties  affected  by  the  temperature  of  the  medium  and 
the  loss  of  coolant  from  the  system? 

(b)  The  shock  strength  and,  hence,  the  acceleration  of 
the  wall  Is  unknown;  how  does  the  Inertial  constraint  Influence 
the  deformation  volume  of  the  reactor  vessel? 

(c)  The  dynamic  tensile  strength  of  the  vessel  wall 
varies  significantly  with  gross  plastic  strain  and  the  rate 
of  strain.  The  shock  loading  and  Inertial  constraints  are 
unknown;  what  are  the  strain  and  strain-rate  functions? 

(d)  The  explosion  gases  mix  with  the  metal-coolant-gas 
medium  surrounding  the  core.  The  character  of  the  Interface 
between  the  explosion  fireball  and  the  surrounding  medium  Is 
unknown;  what  are  the  millisecond  chemical  and  thermodynamic 
mechanisms  that  govern  the  heat-sink  capacity  of  the  medium? 

(e)  The  TMT  Internal  blast  pressure  and  propellant- type 
gas  pressures  depend  upon  the  heat-sink  capacity  of  the 
reactor  Internals  and  the  volime  Into  which  the  equlllbrlim 
gases  expand.  But  the  heat-sink  and  gas-expansion  quantities 
are  unknown  functions  of  unknown  functions,  etc.,  with  coolant 
content  appearing  as  a  significant  variable  In  all. 
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5.3  State  of  the  Art.  Upon  oonalderatlon  of  the  above 
questions  and  some  reflection.  It  was  not  difficult  to 
perceive  the  extremely  complex  functional  form  of  the  pseudo 
parameters  that  characterize  containment.  The  governing 
chemical,  thermodynamic,  and  physical  phenomena  were,  for 
the  most  part,  micro  or  milliseconds  mechanisms  that  did 
not  lie  In  the  realm  of  classical  knowledge.  In  short,  when 
the  program  was  started,  the  state  of  the  art  would  not  permit 
a  quantitative  characterization  of  the  pseudo-parameters  via 
available  theory  and  empiricism.  This  meant  that  one  and  only 
one  method  of  solution  could  be  employed:  Salnt-Venant 's 
semi-inverse  method  of  classical  mechanics. 

5.4  Salnt-Venant 's  Semi-Inverse  Method.  Since  teleological 
solutions  (those  directed  to  a  specific  purpose)  do  not  require 
that  everything  be  known  about  all  of  the  parameters,  bound 
and/or  lumped  parameter  solutions  frequently  provide  Important 
general  Information  In  a  restricted  domain  of  Interest.  One 

of  the  most  useful  algorisms  for  obtaining  such  solutions  Is 
Salnt-Venant 's  semi-inverse  method  of  classical  mechanics. 

Here,  upon  Identifying  and  characterizing  the  principal  param¬ 
eters  to  the  extent  In  which  they  are  known  (both  empirically 
and  theoretically),  a  synthesis  of  these  parameters  Is  used  to 
construct  an  Incomplete  functional  form  of  the  desired  solution. 

This  embryonic  functional  form  constitutes  certain  necessary 
conditions  that  the  complete  solution  must  satisfy.  In  general, 
however,  the  necessary  conditions  expressed  In  the  embryonic 
function  are  not  sufficient  to  guarantee  the  complete  solution; 
this  must  be  achieved  through  the  Judicious  Imposition  of 
boundary  conditions.  Depending  upon  the  nature  of  the  problem, 
the  boundary  conditions  may  be  determined  from  further  theoreti¬ 
cal  considerations,  experimental  data,  or,  very  often,  a  combina¬ 
tion  of  the  two. 
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Of  cours#.  It  may  be  correctly  Inferred  that  the  eeml- 
Inverse  method  will  not  yield  a  solution  If  the  embryonic  form 
cannot  be  adequately  determined.  If  the  character  of  the 
multi-parameter  problem  Is  extremely  complex,  the  state  of  the 
art  will  frequently  not  yield  a  solution  without  numerous 
experimental  data.  This  was  the  case  with  containment:  It 
meant  that  the  necessary  boundary  conditions  could  be  established 
only  through  fundamental  and  extensive  experiments  In  scaled 
models  of  reactor  vessels. 

5.5  Hopklnson  Scaling.  There  exist  many  schemes  for  predicting 
the  behavior  of  a  prototype  from  experiments  In  scaled  models; 
see,  e.g.,  H.  0.  Snay's  specific  treatment  relating  to  underwater 
explosion  phenomena  In  reference  (m)  and  H.  L.  Langhaar's  moire 
general  treatment  In  reference  (n) .  A  study  of  the  various 
methods  showed  Hopklnson  scaling  (also  known  as  Iso-veloclty  or 
cube-root  scaling)  to  be  the  most  appropriate  for  the  containment 
problem. 

Briefly,  In  Hopklnson  scaling,  linear  dimensions  of  the 
prototype  and  model  are  related  by  a  constant  scale  factor  X, 
but  velocity  and  the  acceleration  due  to  gravity  do  not  scale. 
This  fundamental  choice  determines  the  character  of  all  other 
Hopklnson  quantities.  Classical  model/prototype  ratios  of 
quantities  basic  to  containment  are  listed  below: 

Velocity,  pressure,  stress,  strain  .  unity 

length,  time,  displacement,  specific  Impulse  ...  X 

area,  surface  force.  Inertia  .  X^ 

mass,  weight,  vol\ime,  momentum.  Impulse .  X^ 

acceleration,  time  rate  of  strain .  X”^ 

A  fundamental  requirement  of  all  scaling  methods  Is  that 
phenomena  of  Interest  In  the  prototype  be  predictable  from 
experiments  conducted  In  the  model,  and  such  will  be  the  case 
If  the  scaling  character  of  each  phenomenon  Is  known.  Upon 
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relating  this  requirement  to  the  containment  problem,  the 
difficulty  in  obtaining  a  solution  became  more  evident:  for 
example,  not  only  was  the  quantitative  chemical  and  thermo¬ 
dynamic  character  of  the  close-in  TNT-coolant  reaction  generally 
unknown,  its  scaling  properties  were  unknown  as  well.  Specifi¬ 
cally,  this  scaling  question  was  whether  the  TNT-coolant 
mechanism  in  the  prototype  and  model  would  vary  time  wise  with 
X  as  in  Hopkins on  scaling,  or  be  time  dependent  in  another  way. 
More  generally,  there  existed  the  question:  What  are  the 
scaling  properties  of  excursion  loading,  considering  that  the 
reactor  vessel  deforms  dynamically  and  gross-plastlcally? 

5.6  Conceptual  Design  of  Models.  Power  reactor  vessels  and 
many  others  are  complex  welded  structures  that  are  basically 
cylindrical  in  shape  with  length  ^  two  diameters.  The  end 
closures  are  of  various  types,  and  the  vessel  wall,  including 
the  closures,  is  Interrupted  by  numerous  nozzles  to  accommodate 
pipes,  control  devices,  etc.  Aside  from  housing  the  core, 
blanket,  and  control  mechanisms,  the  reactor  vessel  is  filled 
or  partially  filled  with  a  circulating  coolant.  Now  as  dis¬ 
cussed  previously,  the  spectrum  of  excursion  accidents  is  quite 
large  in  both  magnitude  and  flux  of  the  energy  released :  it 
ranges  from  small  accidents  to  thousands  of  megawatt-seconds 
occurring  in  times  from  a  few  microseconds  to  seconds  or  more. 
Furthermore,  the  accidents  of  this  spectrum  can  take  place 
ansrwhere  in  the  coolant  loss  domain  from  no  loss  to  a  complete 
loss . 

The  range  of  excursions,  the  conditions  under  which  they 
may  occur,  and  the  multiplicity  and  complexity  of  the  principal 
parameters  made  it  clear  that  the  general  containment  problem 
was  one  of  Inordinate  scope  and  difficulty.  Rather  than  attempt 
an  immediate  solution  to  the  general  problem,  it  seemed  prudent 
to  seek  first  a  basic  upper-bound  solution  -  a  solution  for 
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containing  a  postulated  excursion  vinder  the  most  adverse 
conditions  In  a  vessel  of  elementary  shape,  A  priori,  the 
"most  adverse  conditions”  for  an  excursion  of  given  magnitude 
are: 

(a)  that  the  excursion  occur  at  the  greatest  energy  flux, 
l.e,,  that  the  energy  be  released  In  the  shortest  time  and 

(b)  that  the  reactor  vessel  be  filled  with  liquid  coolant 
or  suffer  no  loss  of  coolant. 

For  obvious  reasons  the  right  circular  cylinder  was  chosen 
as  the  most  convenient  and  useful  shape  for  the  reactor  vessel 
models.  Early  experiments  with  TNT  In  water  and  air  showed  that 
the  near-rupture  deformation  pattern  (see  Figure  6.1)  Is  always 
restricted  to  about  two  diameters;  hence  the  majority  of  the 
models  possessed  length  equal  to  two  dianeters.  Since  the 
density  and  compressibility  of  water  are  only  slightly  different 
from  those  of  other  known  liquid  reactor  coolants  and  since. 

In  the  containment  domain,  TNT  shock  attenuation  Is  believed  to 
be  negligibly  sensitive  to  these  differences,  water  was  employed 
as  the  model  reactor  coolant. 

The  basic  upper-bound  containment  solution,  then,  was  for 
the  case  of  the  TNT  explosion  occurring  In  a  right  cylindrical 
vessel  with  length  equal  to  two  diameters  and  filled  with  liquid 
coolant.  Once  the  basic  upper-bound  solution  was  In  hand,  the 
rationale  was  that  It  could  be  modified  to  provide  safe  and 
reasonable  solutions  for  less  adverse  accident  conditions  and 
for  vessels  fitted  with  weldments,  nozzles,  and  other  dis¬ 
continuities  . 

5.7  Summary.  The  principal  parameters  of  containment,  both 
classical  and  pseudo,  were  set  forth,  and  some  fundamental 
questions  were  raised  to  Indicate  the  difficulty  of  quantitative 
parametric  characterization.  An  upper-bovmd,  Ivmiped-parameter 
solution  via  Saint -Tenant 's  semi-inverse  method  was  postulated 
as  the  only  path  that  would  yield  the  desired  Containment  Laws. 
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The  senl-lnyerse  method,  Hopklnson  scaling,  and  conceptual 
models  were  discussed,  and  case  conditions  for  the  basic  upper- 
bound  containment  solution  were  specified. 

The  semi-inverse  method  Is  pursued  In  subsequent  chapters 
via  the  following  steps: 

(a)  formulation  of  the  semi-inverse  embryonic  solution 
(l.e.,  an  embryonic  Containment  Law  format) 

(h)  conduct  of  an  extensive  experimental  program  employing 
explosives  and  propellants  to  simulate  nuclear  exctirslons  In 
Hopklnson  models  of  reactor  vessels 

(c)  reduction  of  experimental  data  to  satisfy  boundary 
condition  requirements  of  the  embryonic  solution 

(d)  modification  of  the  basic  upper-boxmd  solution  to 
provide  safe  and  reasonable  solutions  for  the  less  adverse 
conditions  of  the  accident  spectrum  and  for  vessels  fitted  with 
weldments,  nozzles,  and  other  discontinuities. 
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CHAPTER  6 
EMBRYONIC  FORMAT 

The  flrat  step  In  the  Semi-Inverse  Method  of  solution 
Is  to  synthesize  everything  that  Is  known  (or  seems  reasonable) 
about  the  problem  and  to  oonstruot  an  embryonic  or  trial  solu¬ 
tion  from  the  synthesis. 

6,1  Preliminary  Experiments.  An  examination  of  available  con¬ 
tainment  Information  showed  that  even  an  adequate  trial  format 
of  a  basic  containment  law  could  not  be  written  down  without 
some  preliminary  experiments.  Among  the  first  of  these  were 
tests  employing  pentollte  In  right-cylindrical  model  vessels 
filled  with  water  or  air  at  atmospheric  pressure.  The  principal 
results  were  the  "localization  factors"  shown  In  Figure  6,1. 
(Figures  and  tables  appear  at  the  end  of  each  chapter.)  Here 
It  Is  seen  that  the  character  of  the  shock  Is  such  that  near¬ 
rupture  deformation  Is  localized  to  about  two  radii  for  air  and 
four  radii  for  water.  Hence,  for  vessels  with  Inside  length  L 
and  Inside  radius  Rj[  such  that  L  >  4Rj^,  it  was  Indicated  that 
the  explosion  containment  potential  Is  not  an  appreciable  function 
of  vessel  length.  A  priori,  this  w<mld  be  true  Irrespective  of 
the  loss  of  reactor  coolant  In  any  degree. 

In  other  preliminary  experiments  employing  pentollte  In 
water-filled  models,  dilations  were  photographed  with  A  Imeloian- 
Whltley  Dynafax  camera  operating  at  23,000  frames  per  seeMid. 
Typical  vessel  data,  test  conditions,  and  sequential  frames  are 
shown  In  Figure  6.2.  The  total  dilation  time  was  about  600 
microseconds;  at  40  microseconds  per  frame,  the  motion  was 
recorded  In  13  frames.  Radial  deformation  and  strain  of  the 
external  wall  at  the  mld-merldlan  Is  plotted  versus  time  In 
Figure  6.3.  From  these  data,  radial  velocity  and  strain  rate 
versus  time  are  plotted  In  Flgvire  6.4;  radial  velocity  and 
strain  rate  versus  radial  strain  are  plotted  In  Figure  6.3. 
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The  geometries  of  the  curves  In  Figures  6.3,  6.4,  and 
6.5  are  singularly  significant.  The  def ormatlon-tlme  plot 
can  be  Interpreted  as  parabolic.  If  this  Is  done.  It  follows 
that  the  velocity-strain  plot  Is  parabolic  and  that  the 
velocity-time  plot  Is  linear.  Prom  D'Alembert's  Principle, 
this  means  that  the  force  acting  to  bring  the  mld-merldlan 
wall  to  rest  was  constant  throughout  the  entire  deformation. 
Several  Important  consequences  of  these  fladlngs  were: 

(a)  Specifically,  the  data  provided  Important  parametric 
Information  for  the  embryonic  containment  law 

(b)  More  generally,  new  Insight  was  provided  Into  the 
response  of  pressure  vessels  to  explosion  loading  In  the  con¬ 
tainment  domain  (the  response  domain  In  which  containment  occurs), 
and 

(c)  Mew  Insight  and  new  corroborative  data  were  added  to 
the  meager  Information  available  on  the  effects  of  very  rapid 
strain  rates  upon  the  strength  of  ductile  steels.  These  data 
were  obtained  for  a  mld-merldlan  ring  that  deformed  gross - 
plastically  under  the  action  of  an  apparently  constant  retarding 
force.  Only  this  force,  a  function  of  the  dynamic  strength  of 
the  material,  and  the  D'Alembert  term  determined  the  motion. 

Hence,  these  rate -of -strain  data  are  Inherently  free  of  the  end. 
Inertial,  and  short  loading  time  effects  so  common  to  dynamic- 
tensile  and  Impact  experiments. 

The  results  of  the  preliminary  experiments  together  with 
available  explosions  phenomenology  motivated  and  permitted 
formulation  of  the  embryonic  solution,  first  for  the  general 
case  and  then  for  a  specific  case. 

6.2  General  Case.  A  closed  vessel  Is  subjected  to  Internal 
explosion  loading.  The  vessel  wall  Is  a  strain-hardening  steel 
that  flows  gross-plastlcally.  It  Is  assumed  that  the  wall 
material  Is  of  uniform  constant  mass  density  and  that  the  wall 
remains  a  continuum. 
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We  frame  the  vessel  in  an  orthogonal  coordinate  system 
for  which  the  Lagranglan  spatial  coordinates 

^lo  >  1''1>2,3 


denote  the  original  position  (l.e.,  the  position  at  event  time 
t  =  O)  of  all  points  In  the  vessel  wall  Including  the  surface. 
Since  the  Eulerlan  coordinates 


xi  *=  xi(xio,  t) 


are  at  least  of  class 


and  the  Jacobian 


d  (x^,  x«,  x^) 

- 7 — ^ ^ ^ — r"  •=  J 

h  ^30^ 

Is  everywhere  positive  in  the  wall,  the  EuDerlsn  and  Lagranglan 
Integrals  are  related  In  the  usual  way  by 


f  f  (xi)  dxi  /  f  (xj )  J  dx^, 


where  the  dxj^  and  dXj^^  are  products,  rcspecci\rely,  of  the  dxi 
and  dx^Q. 
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Upon  detonation  of  the  explosive,  the  shock  wave  propa¬ 
gates  through  the  media  surrounding  the  charge  and  strikes 
the  Internal  wall;  whereupon  the  wall.  It  Is  hypothesized  here, 
responds  with  a  virtually  Instantaneous  Initial  velocity 

The  Instantaneous -velocity  hypothesis  Is  reasonable  since  the 
effective  shock  loading  time  Is  an  order  of  magnitude  less  than 
the  gross  deformation  time.  Velocity-time  and  velocity-strain 
Integrations  of  the  second  law  of  motion,  respectively,  yield 


P  (Mo,  t)  •  dA 


lo)  ^^10 


w 


h  dt 


i/ 

tr 


Vo(='lo)'‘='lo 


)  <1*10 


(6.1) 


(6.2) 


where  the  symbology  Is 


▼o 


t 

P 


W 


g 

A,  V 
D 


Initial  velocity  of  vessel  wall,  ft/sec 
duration  of  shook  Interval,  sec 
unit  vector  In  direction  of  dA 
vector  shock  pressure,  psfg 
wall  mass  density,  Ib/ft^ 

p 

acceleration  due  to  gravity,  ft/sec 
wall  Internal  area  and  volume,  ft^,  ft^ 
wall  deformation  energy  density,  ft-lb/ft^ 
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Equation  (6.1)  states  that  the  total  Impulse  delivered  to  the 
vessel  wall  Is  equal  to  the  total  Initial  monentiim  of  the  wall; 
(6.2)  states  that  the  total  work  done  In  deforming  the  wall  It 
equal  to  the  total  Initial  kinetic  energy  of  the  wall.  Neither 
the  momentiun  nor  the  kinetic  energy  Is  conserved,  but  the  con¬ 
servation  of  mechanical  energy  Inherent  In  (6.2)  Is  valid. 

6.3  Specific  Case.  While  the  general  relations  expressed  In 
equations  (6.1)  and  (6.2)  possess  a  conceptually  desirable  form, 
they  will  not  yield  a  quantitative  containment  solution  In  terms 
of  currently  available  Information.  We  turn  then  to  a  less 
general  solution:  one  that  Is  further  restricted  to  a  water- 
fllled  rlght-cyllndrloal  vessel  for  which 

(a)  the  Initial  values  of  the  Internal  and  external  radii, 
Rj^  and  R^,  the  Internal  length  L,  and  the  wall  thickness  h^  are 
related  by 


L  i  10  s  RiAo  ^  ^0  ^o  “  ^e  "  ^1 

(b)  Rj^  and  the  explosive-charge  radius  R^  are  related  by 


Rl  ^  3Rc 


(c)  the  water  Is  In  a  liquid  state  attendant  normal 
operating  conditions* 


N  It  should  be  noted  that  the  analysis  and  experimental  work 
to  follow  assumes  the  water  to  be  at  room  temperature  and 
atmospheric  pressure  rather  than  at  elevated  values.  This, 
however,  does  not  constitute  a  significant  restriction.  The 
elevated  temperatures  and  pressures  normally  associated  with 
sodium  or  water  reactors  should  not  affect  appreciably  the 
shock  Impulse  delivered  to  the  vessel  wall.  Any  elevated 
pressure  would  alter  only  the  elastic  strain  energy  of  the 
vessel,  but  this  Is  considered  negligible  In  the  gross-plastic 
deformation  domain. 
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(d)  the  vessel  wall  is  hcxnogeneous  and  isotropic 

(e)  no  constraints  other  than  the  tensile  and  shearing 
constraints  in  the  wall  act  to  retard  gross  plastic  deforma¬ 
tion  of  the  vessel 

(f )  the  explosive  charge  is  compact  (length  «  diameter) 
and  is  detonated  at  the  centroid  of  the  vessel 

(g)  there  are  no  nozzles,  welds,  or  other  wall  discon¬ 
tinuities  . 

The  above  restrictions  contain  three  Inequalities.  The 
first  of  these,  L  ^  4Rj^,  is  motivated  by  the  "localization 
factor"  for  water  illustrated  in  Figure  6.1.  It  is  mentioned 
again  that  the  character  of  the  shock  is  such  that  deformation 
is  localized  to  not  more  than  four  radii. for  either  water  or 
air  surrounding  the  charge.  Hence,  for  L  ^  4Rj^,  the  explosion 
containment  potential  of  the  reactor  vessel  simulant  is  not  an 
appreciable  function  of  its  length.  The  lower  limit  of  the 
second  inequality,  Ri/ho  ^  10,  is  an  accepted  value  for  thin- 
wall  theory;  the  upper  limit,  Ri/hQ  ^  40,  is  the  maximum  value 
selected  for  this  study.  The  third  inequality,  Rj  ^  3Rc>  Is 
invoked  because  of  the  extent  to  which  hydrodynamic  phenomenology 
for  the  subject  problem  is  unknown  closer  to  the  charge.  It 
is  Important  to  note  that  since  specifications  of  the  large 
majority  of  reactors  can  be  interpreted  to  satisfy  the  above 
restrictions  on  configuration,  no  loss  in  generality  accrues 
in  this  respect. 

Consistent  with  the  cylindrical  coordinate  system,  we  put 

^10  ”  ’^o^  ^0  ”  ^30  ~ 

and  if  we  place  the  system  origin  at  the  centroid  of  the 
cylinder  with  major  and  z  axes  colinear,  it  follows  from 
symmetry  that  equations  (6.1)  and  (6.2)  reduce,  respectively,  to 
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(6.3) 


(6.4) 


where  r  is  the  unit  vector  In  the  radial  direction. 


Vu  =  TT  (r|  -  Ri)  (6.5) 

is  the  volume  per  unit  length  of  cylinder,  and 

®h  “ 

is  the  average  value  of  D(Zq)  across  the  wall  thickness  h, 

6.4  Unit-Length  Ring.  We  elect  at  this  point  to  consider 
the  response  of  a  single  unit-length  ring  located  at  the  mid- 
merldlan  of  the  cylinder  in  the  (r,  0,  o)  plane  wUbh  rauJli 
Rj[^  and  Rg.  Upon  detonation  of  the  explosive  charge,  this  ring 
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being  closest  to  the  charge  centroid,  receives  the  greatest 
Impulse,  undergoes  the  largest  deformation,  and  Is  subjected 
to  the  highest  stress.  In  fact,  the  hoop  stress  In  the  mld- 
merldlan  ring  constitutes  an  excellent  containment  Index,  since 
experimental  evidence  Indicates,  as  shown  In  Plgvire  6,6,  that 
marginal  ruptxire  for  ductile  cylinders  always  results  from  the 
propagation  of  a  fissure  that  originates  at  the  mld-raerldlan. 
Invoking  this  rationale,  we  find  that,  for  the  mid -meridian 
ring  of  length  unity,  the  variable  Zq  becomes  the  parameter 
Zq  -  0,  and  equations  (6.3)  and  (6.4)  take  the  forms 


2tt  Ri 


dt 


g  ^u  ^o  ( o) 


(6.6) 


Dh(0)  ■  ^  Vo^(o) 


(6.7) 


The  simultaneous  solution  of  equations  (6.5)»  (6.6),  and  (6.7) 
yields 


p(0,t)  dt 


y 


2g  Rl 


I>h(o) 


(6.8) 


p 

where  I  Is  the  specific  Impulse  (lb  sec/ft  ).  Thus  the  functions 
of  Interest  at  this  stage  of  the  embryonic  solution  are  I,  the 
specific  Impulse,  and  Dh,  the  deformation  energy  density. 
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6.5  Specific  Impulse.  An  eetabllahed  power  law  relating 
charge  weight  W  and  charge  distance  Rj[  to  specific  Impulse 
for  explosions  In  free  water  Is  given  by  Cole  In  reference  (l) 
as 


If  » 


(6.9) 


where  the  nomenclature  Is 


If 

W 

Ri 

K,  $ 


p 

specific  Incident  Impulse  In  free  water,  lb  sec/ln 
explosive  charge  weight,  lb 
distance  from  charge  center,  ft 
constants 


The  quantities  K  and  0  are  treated  as  constants  In  the  fx*ee- 
water  domain  for  a  given  explosive.  For  pentollte,  0  «  1.05 
In  the  broad  spectmm  of  conventional  undeznrater  applications; 
for  the  narrow  containment  domain,  however.  It  was  found  that 
the  classical  acoustic  law  (0  «  l)  can  be  assumed  with  no  loss 
In  generality.  Thus  equation  (6.9)  reduces  to 


If 


(6.10) 


Of  course, the  total  Impulse  of  (6.10)  will  not  be  absorbed  by 
the  vessel  wall;  a  portion  will  be  lost  due  to  reflection  and 
the  complex  phenomena  attendant  to  all  such  highly  transient 
high-pressure  Interfaces.  To  account  for  these  losses,  as  well 
as  a  difference  In  units,  we  hypothesise  that  I  (lb  sec/ft^), 
the  Impulse  absorbed  by  the  vessel  wall.  Is  related  to  the  free- 

p 

water  Impulse  Ij.  (lb  sec/ln  )  by  an  "efficiency  factor"  function 
^  such  that 
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I  =  144  If/? 


(6.11) 


It  seemB  reasonable  to  expect  that  ?  Is  a  function  of  the  wall 
distance  from  the  charge  and  the  D'Alembert  and  tensile  con¬ 
straints  In  the  wall.  More  specifically,  we  know  that  the 
classical  hoop  stress  Is  a  function  of  Rj^/h©.  Consistent  with 
the  semi-inverse  method  and  the  teleological  character  of  the 
desired  solution,  we  Incorporate  this  rationale  Into  the  trial 
format  and  write 


?  -  ?  (RlAo) 


(6.12) 


whereupon  the  Impulse  absorbed  by  the  vessel  wall  becomes 


I  -  l44  KW^/3^1? 


(6.13) 


6.6  Deformation  Energy.  The  true  stress-classical  strain 
diagram  Is  quite  linear  for  the  case  of  a  ductile  steel 
deformed  statically  In  the  gross  plastic  domain.  With 
sufficient  accuracy  for  our  semi-inverse  rationale,  we  can 
write  that 


e 


(6.14) 


36 


NOLTR  63-140 


where  the  nomenclature  is 


ot 


true  stress,  psl 
conventional  yield  stress,  psl 
conventional  ultimate  stress,  psl 
conventional  Instantaneous  strain,  In/ln 
conventional  ultimate  strain  (elongation),  In/ln 


For  the  case  of  gross  plastic  deformation  occurring  at  rapid 
rates  of  loading,  the  stress  a  in  a  strain-hardening  steel  like 
304  stainless  Is  at  least  a  function  of  the  strain  c  and  the 
strain  rate  e.  The  stress  may  also  be  a  function  of  V  and 
higher  derivatives  of  e  with  respect  to  time;  but  since  little 
or  nothing  Is  known  about  these  higher  order  effects.  It  will 
be  assumed  here  In  the  usual  way  that  the  stress  function  Is 


a  =  a  (®,  (6.15) 

l.e.,  the  stress  Is  a  function  of  the  strain-hardening  and 
strain-rate  mechanisms  alone.  For  deformations  that  occur  at 
continuously  decreasing  rates  of  strain,  the  strength  of  the 
material  decreases  continuously  with  strain  as  a  function  of 
the  strain-rate  mechanism  and  Increases  continuously  with  strain 
as  a  function  of  the  strain  hardening  mechanism. 

It  Is  Interesting  to  entertain  the  net  Influence  of  the 
strain-hardening  and  strain-rate  mechanisms  Implied  In  Figures 
6.3,  6.4,  and  6.5.  In  reference  (o),  Tardlf  and  Erickson  report 
the  work  of  a  number  of  Investigators  on  the  strain-rate  enhance¬ 
ment  of  the  yield  point  in  mild  steel.  The  data  are  presented 
graphically  with  the  ratio 

R  K  d^amlc  yield  stress 
“  static  yield  stress 
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on  the  ordinate  and  strain  rate  on  the  abscissa.  Among  the 
results  are  those  of  M.  J.  Manjolne  who  fo\md  that  a  strain 
rate  of  1300  per  second  produces  the  ratio  R  w  3.  Manjolne 
also  found  that  the  ratio  R  Is  quite  stable,  varying  exponentially 
from  2-1/2  to  3  for  strain  rates  between  200  and  1300  per  second. 
For  the  stainless  steel  and  loading  conditions  of  Figures  6.2, 
his  results  Imply  an  Increase  In  the  yield  point  from  30,000 
psl  to  90,000  psl,  a  value  that  Is  equal  to  the  true  static 
stress  at  37  per  cent  elongation. 

If  Manjolne's  results  are  accepted  for  the  moment,  this 
means  that  when  the  vessel  started  to  deform  at  a  strain  rate 
of  1300  per  second  and  when  It  came  to  rest  with  a  37  per  cent 
elongation,  the  hoop  stress  at  the  mld-merldlan  was  the  same 
value:  90,000  psl.  This  coincidence  and  the  stability  of  the 
strain-hardening  and  strain-rate  effects  are  consistent  with 
the  finding  In  section  6.1  that  the  force  acting  to  bring  the 
mld-merldlan  wall  to  rest  was  constant  during  the  deformation 
of  this  particular  vessel.  It  Is  Interesting  to  entertain  a 
hypothetical  hoop  stress -strain  diagram  that  might  be  generated 
from  such  loading  conditions.  This  Is  done  In  Figure  6.7,  and 
a  comparison  Is  made  with  statically  generated  diagrams.  From 
a  general  point  of  view,  however,  It  must  be  noted  that  strain 
rate,  for  the  conditions  under  discussion,  varies  Inversely  with 
the  scale  factor,  X.  This  means  that  In  larger  vessels  the 
strain-rate  enhancement  of  tensile  strength  would  be  diminished 
and  that,  of  course,  the  mld-merldlan  force  would  no  longer  be 
a  constant. 

Manjolne  also  presented  strain-rate  data  for  mild  steel  In 
reference  (p).  Bodner  and  Symonds,  In  reference  (q),  fitted  a 
curve  to  these  data  and  found  that 


(6.16) 
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where  the  nomenclature  Is 

0(1  ....  dynamic  yield  stress,- psl 

Oy  ....  static  yield  stress,  psl 
4  ....  strain  rate,  per  second 

Johnson  In  reference  (r)  reported  results  of  explosion 
tests  with  small  unconstrained  rings  of  various  materials. 

The  experimental  techniques  that  he  employed  eliminated 
longitudinal  and  shear  stresses  and  ensured  a  virtually 
Instantaneous  wall  velocity  with  no  further  loading  once 
the  ring  began  gross  motion.  From  D’Alembert’s  Principle 
he  determined  the  time  dynamic  hoop  stress  In  the  ring  as 
a  function  of  strain,  and  for  several  ductile  materials  he 
found  that  this  stress,  aside  from  being  elevated  with  strain 
rate,  was  linear  with  strain.  Of  course,  there  exist  other 
strain-rate  data  that  could  be  listed  here,  but  the  above 
examples  and  discussion  are  sufficient  to  Indicate  the 
rationale  attendant  the  development  of  a  deformation  energy 
function  for  strain  that  Is  dynamic  and  gross  plastic. 

A  priori  and  specifically  from  the  relations  (6,l4)  and 
(6.15)>  we  try  the  embryonic  form 


D^{0)  -  144  at  c  (i^) 


(6.17) 


where  la  the  Initial  strain  rate  corresponding  to  the  Initial 
wall  velocity  v^.  The  simultaneous  solution  of  equations  (6.8), 
(6.13),  and  (6.17)  yields 


W  -  [  ^2  wat*  (Re  '  «l)^ 


(6.18) 
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where 

-  ’l  («o)  .  ’2  -  25s|^  “  »2  (RiAo)  (6-19) 


Equation  (6.18)  conetltutee  the  embryonic  containment  law 
format  where  an  extensive  experimental  program  Is  required 
to  establish  the  deformation  energy  function  and  the 
efficiency  factor  function  Yg, 
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NOTES: 

VESSEL  MATERIAL  AND  SHAPE;  10"  ID,  20"  LENGTH,  1/4"  WALL,  ROLLED 
FROM  304  STAINLESS  STEEL  PLATE,  WELDED  ALONG  LONGITUDINAL  SEAM. 

TEST  CONDITIONS;  VESSEL  FILLED  WITH  WATER  AT  700F  AND  ATMOS¬ 
PHERIC  PRESSURE,  200  GRAMS  PENTOLITE  DETONATED  AT  CENTROID. 

PHOTOGRAPHY:  VESSEL  DILATION  RECORDED  BY  BECKMAN  -  WHITLEY 
DYNAFAX  CAMERA  OPERATING  AT  25,000  FPS.  FRAMES  SHOWN  ABOVE 
TAKEN  AT  40 -MICROSECOND  INTERVALS,  TOTAL  DILATION  TIME -600 
MICROSECONDS  (15  FRAMES). 


FIG.  6.2  EXPLOSIVE  DILATION  OF  MODEL  VESSEL 
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FIG.  6.5  VELOCITY -STRAIN  RATE-STRAIN  PLOT  OF  EXPLOSIVELY  DILATED  VESSEL 
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19  GRAMS  OF  PENTOLITE  IN  WATER 
PRODUCED  THIS  MARGINAL  RUPTURE 


TYPICAL  MARGINALLY  RUPTURED  VESSEL  SHOWS  THAT 
FISSURES  ORIGINATE  AT  THE  MID-MERIDIAN  AS  EXPECTED. 
WALL  MATERIAL  WAS  CENTRIFUGALLY  CAST  304  STAINLESS 
STEEL.  ORIGINAL  DIMENSIONS  WERE  5"  ID,  10"  INTERNAL 
LENGTH,  1/8"  WALL  THICKNESS. 


FIG.  6.6  DEFORMED  VESSEL  SHOWING  ORIGIN  OF  RUPTURE 
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FIG.  6.7  COMPARISON  OF  DYNAMIC  AND  STATIC  STRESS  -  STRAIN  DIAGRAMS 
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CHAPTER  7 

EXPERIMENTAL  PROGRAM  AND  RESULTS 

The  following  experimental  program  provided  necessary 
boundary  conditions  for  evaluating  the  deformation  energy  and 
efficiency  factor  functions  of  the  embryonic  containment  law 
format. 

7.1  Idealized  Models.  The  major  effort  of  the  experimental 
program  was  directed  to  explosively  dilating  model  reactor 
vessels  without  nozzles  or  significantly  restrictive  end  con¬ 
straints.  In  each  of  these  experiments: 

(a)  the  model  vessel  was  provided  with  a  weak  water-tight 
closure  at  the  lower  end  (the  upper  end  was  open) 

(b)  the  vessel  was  filled  with  water  at  atmospheric  pres¬ 
sure  and  70°P 

(c)  a  compact  (length  ~  diameter)  pentollte  charge  was 
detonated  at  the  centroid  of  the  vessel. 

Although  the  models  were  limited  to  convenient  sizes  for 
laboratory  testing,  an  ample  spectrum  of  sizes  and  wall  thick¬ 
nesses  was  Investigated,  as  listed  In  the  following  table: 

ID  (In)  Length  (in)  Wall  Thickness  (in) 


5 

10 

1/16,  1/8,  1/4 

5 

20 

1/8 

10 

20 

1/8,  1/4,  1/2 

15 

30 

1/2 

20 

40 

1/4,  1/2,  1 

20 

6o 

1 

To  provide  containment  law  knowledge  for  a  wide  range  of 
ductile  materials,  model  vessels  were  fabricated  from  1020  mild 
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steel,  212  flange  steel,  several  types  of  stainless  steel,  and 
one  type  of  aluminum.  These  materials  and  their  mechanical 
properties  are  given  In  the  table  below: 


Material 

Yield 

Stress 

(psl) 

Ultimate 

Stress 

(psl) 

Ultimate 

Strain 

(In/ln) 

304  Stainless  Steel 

30, 000 

65,100 

0.60 

304  Stainless  Steel* 

43,000 

83,800 

0.60 

4l0  Stainless  Steel 

66,000 

94,300 

0.26 

4l0  Stainless  Steel 

151,500 

197,900 

0.08 

212  Flange  Steel* 

46,500 

76,200 

0.27 

1020  Mild  Steel* 

42,000 

66,900 

0.30 

6o61-T6  Aluminum 

40, 000 

45,000 

0.15 

The  stress  and  strain  values  are  averages  obtained  from  nimierous 
tensile  test  specimens  of  each  material.  Vessels  made  of 
materials  denoted  by  an  asterisk  (*)  were  constructed  of  rolled 
plate  with  one  or  two  longitudinal  welded  seams.  In  general, 
the  light-wall  vessels  were  made  with  a  single  weld;  the  heavy- 
wall  vessels  were  made  with  two  welds.  Welds  represented  a 
departiire  from  the  Idealized  vessel  concept,  but  special  pre¬ 
cautions  were  taken  to  ensure  that  the  welds  were  capable  of 
withstanding  the  same  stress  and  deformation  as  the  parent  metal; 
thus  the  welded  vessels  could  be  considered  as  Idealized  models. 
The  precautions  were: 

(a)  the  welds  were  continuous  and  full  penetration 

(b)  the  weld  material  possessed  essentially  the  same 
mechanical  properties  as  the  parent  metal 

(c)  the  vessels  were  annealed  to  relieve  any  residual 
and  thermal  stresses  resulting  from  the  rolling  and  welding 
processes. 
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By  virtue  of  these  precautions,  welds  of  equal  or  perhaps 
better  quality  than  those  expected  in  actual  reactor  vessels 
were  ensured.  The  remaining  stainless  steels  in  the  above 
table  were  in  the  form  of  centrifugally  cast  tubing  (seamless), 
and  the  sole  aluminum  was  in  the  form  of  extruded  seamless 
tubing.  Figure  7.1  shows  typical  test  vessels  of  the  seamless 
and  welded  varieties. 

The  test  procedure  followed  throughout  the  program  con¬ 
sisted  of  explosively  dilating  a  minimum  of  four  vessels  for 
each  given  material  and  size.  In  this  manner  it  was  possible 
to  approach  and  determine  the  condition  of  marginal  containment. 
(Marginal  containment  is  consistent  with  the  maximum  charge 
weight  that  can  be  detonated  in  a  vessel  without  causing  rupture.) 
In  terms  of  experimental  procedure,  three  vessels  were  employed 
to  yield  intermediate  values  of  terminal  radial  strain  (the  at- 
rest  mid -meridian  strain)  correlative  with  charge  weight.  The 
fourth  yielded  the  maximum  radial  strain  correlative  with 
marginal  containment,  and  a  25, 000-frame/second  Beckman-Whltley 
Dynafaix  camera  recorded  the  deformation-time  history  during 
dilation.  Treatment  of  these  photographic  data  in  the  manner 
described  in  section  6.1  resulted  in  the  determination  of 
initial  strain  rates  for  the  various  vessels. 

Table  7.1  gives  the  following  data  for  104  experiments  in 
water-filled,  idealized  model  vessels. 

(a)  charge  weight 

(b)  tennlnal  radial  strain  at  the  mid -meridian 

(c)  mode  of  vessel  containment  or  failure 

(d)  initial  strain  rate  as  determined  from  the  high-speed 
photographs . 

-  Figure  7.2  shows  two  10"  ID,  20"  length,  1/4"  wall,  304 
stainless  steel  cylinders  that  typify  deformation  patterns  of 
marginal  containment.  Several  observations  can  be  made  from 
the  data  of  Table  7.1;  they  are; 
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(a)  For  a  given  veesel  size  and  material,  a  high  degree 
of  uniformity  in  the  terminal  radial  strain  data  is  noted  as 
charge  weight  increases  toward  the  marginal  containment  weight. 

This  is  demonstrated  in  Figure  7,3,  which  shows  four  304  stainless- 
steel  cylinders  in  stages  of  increasing  deformation. 

(b)  In  nearly  all  cases  of  marginal  containment,  the 
maximum  radial  strain  was  slightly  less  than  the  correlative 
average  ultimate  strain  given  in  the  materials  table  earlier 
in  this  section.  The  two  exceptions  were  the 

10"  ID,  20"  length,  1/4"  wall,  4l0  annealed  stainless 
steel  vessels,  and  the 

20"  ID,  40"  length,  1/2"  wall,  212  flange  steel  vessels. 

Not  only  did  the  strains  of  these  materials  exceed  the  listed 
average  values,  they  exceeded  the  maximum  ultimate  strains  given 
by  any  of  the  respective  certified  tensile  test  specimens. 

7.2  Effects  of  Welds.  It  is  recalled  from  section  7.1  that 
fabrication  of  the  rolled-plate  vessels  was  such  as  to  ensure 
welds  of  containment  potential  equivalent  to  that  of  the  parent 
material.  The  success  of  this  effort  was  evidenced  in  vessel 
responses  such  as  those  shown  in  Figures  7.2  and  7.3.  It  is 
noted  that  the  deformation  capacity  of  the  welded -plate  vessels 
was  not  less  than  that  of  their  seamless  counterparts.  Although 
the  majority  of  the  welded  vessels  demonstrated  excellent  deforma¬ 
tion  properties,  there  were  cases  of  premature  rupt\ire  due  to 
weld  failures  either  in  the  weld  material  itself  or  in  the  heat- 
affected  zone  adjacent  to  the  weld.  However,  it  is  seen  from 
Table  7.1  that  in  all  cases  the  vessels  were  capable  of  deforming 
without  rupture  to  strains  greater  than  1/3  and,  with  only 
several  exceptions,  to  strains  greater  than  1/2  c^.  From  these 
data  it  was  postulated  that  the  maximum  strain  of  welded  reactor 
vessels  should  be  restricted  to 
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C  S  Cu/3 


7.3  Effects  of  End  Closures.  Of  the  104  experiments  reported 
In  Table  7.1,  nine  were  conducted  in  vessels  fitted  with  rigid 
radial  end  constraints  of  the  type  shown  in  Figure  7.4.  These 
vessels  were  5"  ID,  10"  length,  1/8"  and  1/4"  wall,  304  stainless 
steel  (cast)  cylinders  and  are  denoted  in  Table  7.1  by  ♦♦♦.  It 
is  seen  immediately  from  the  experimental  results  that  the  end 
closures  decreased  substantially  the  strain  at  which  rupture 
occurred:  from  approximately  0.50  to  0.38  In/in.  Prom  these 
data  it  was  postulated  that  the  permissible  strain  for  a  vessel 
with  rigid  end  closures  should  be  restricted  to 

€  ^  6^/2 

It  can  be  expected  that,  for  the  vessels  fittea  with  hemi¬ 
spherical  or  dished  end  closures,  the  permissible  strain  would 
exceed  that  for  the  case  of  rigid  ends. 

7.4  Effects  of  Nozzles.  To  index  the  limiting  effects  of 
nozzles  on  containment,  experiments  were  conducted  under  the 
following  four  conditions: 

(a)  vessel  fitted  with  single  circular  aperture  at  mid¬ 
meridian;  aperture  diameter  1/5  that  of  vessel  diameter 

(b)  vessel  fitted  with  single  closed-end  tube  inserted 
at  mid-meridian  in  aperture  of  same  diameter  as  in  (a);  tube 
secured  to  vessel  wall  by  continuous  fillet  welds  around  cir- 
cvunference  of  insert 

(c)  vessel  same  as  in  (a)  except  that  three  apertures 
were  equally  spaced  around  circumference  at  mid-meridian 

(d)  vessel  same  as  in  (c)  except  that  three  closed-end 
tubes  were  inserted  in  apertures  and  secured  to  vessel  wall  by 
continuous  fillet  welds  around  circumference  of  inserts. 
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Types  (a)  and  (b)  are  shown  undeformed  in  Figure  7.5;  post- 
test  deformations  are  shown  in  Figure  7.6.  Experimental  results 
from  12  tests  are  given  in  Table  7.2.  Althovigh  these  experi¬ 
ments  do  not  represent  a  comprehensive  study  of  nozzles,  they 
permitted  the  following  upper-bound  observations. 

(a)  For  a  single  mid -meridian  aperture  of  diameter  ~  1/5 
the  diameter  of  the  vessel,  the  maximum  permissible  strain  for 
304  stainless  steel  was  about  0.10  In/ln  or  1/6  c^. 

(b)  For  a  nozzle  welded  on  both  sides  of  the  vessel  wall, 
the  maximum  permissible  strain  for  304  stainless  steel  was  0.33 
in/in  or  about  1/2 

(c)  The  presence  of  multiple  apertures  or  nozzles  at  the 
mid-meridian  did  not  affect  appreciably  the  respective  contain¬ 
able  charge  weights  for  the  cases  of  one  aperture  or  one  nozzle, 
althovigh  the  average  clrcvimferentlal  strain  was  less.  However, 
the  local  deformation  and  thickness  of  the  vessel  wall  at  points 
equidistant  from  the  discontinuities  appeared  to  have  the  same 
characteristics  as  the  uninhibited  deformation  of  idealized 
models  for  equivalent  charge  weights. 

On  the  basis  of  these  experiments,  it  was  postulated  that 
the  maximum  permissible  strain  for  vessels  with  nozzles  (no 
constraints  extemal  to  vessel)  should  be 

e  i  eu/2 

There  is  additional  reason  to  believe  that  this  strain  restriction 
is  conservative;  the  integrity  of  the  welded  nozzles  of  the 
models  is  not  considered  as  sound  as  the  flanged  or  forged  nozzles 
found  in  typical  reactor  vessel  construction.  Appendix  B  gives 
additional  data  concemlng  the  strain  pattern  around  an  aperture 
and  nozzle  for  the  two  vessels  shown  in  Figvire  7.6.  These  data 
should  be  useful  to  those  concerned  with  the  detail  design  of 
nozzle  connections  for  reactor  vessels  and  ancillary  piping. 


53 


NOLTR  63-140 


7.5  Loss -of -Coolant  Experiments.  To  determine  the  enhance¬ 

ment  in  explosion  containment  potential  resulting  from  the 
complete  loss  of  liquid  coolant,  a  limited  number  of  experiments 
were  conducted  in  5"  10”  length,  l/l6"  and  1/8"  wall,  304 

stainless  steel  (cast)  vessels  of  the  type  shown  in  Figure  7.4. 
In  each  test  the  vessel  was  filled  with  air  at  atmospheric 
pressure  and  ambient  temperature,  and  the  explosive  charge  was 
detonated  at  the  centroid  of  the  vessel.  Charge  weight  - 
terminal  strain  results  are  given  in  Table  7.3,  and  the  uniform 
Increase  in  deformation  with  charge  weight  for  1/8"  wall  models 
is  shown  in  Figure  7.7.  The  localized  deformation  pattern  and 

a  comparison  with  that  of  water-filled  vessels  were  shown  pre¬ 
viously  in  Figure  6.1. 

7.6  Applicability  of  Free-Water  Explosion  Relations.  In  the 
development  of  the  embryonic  format  (Chapter  6),  it  was  assumed 
that  conventional  free-water  explosion  relations  could  be 
employed  in  characterizing  the  wall  loading.  However,  Cole 

in  reference  (l)  states  that  the  use  of  these  relations  is 
questionable  for  distances  closer  to  the  charge  than  7  to  10 
charge  radii.  In  some  of  the  water-filled  vessel  experiments, 
the  distance  to  the  vessel  wall  was  as  small  as  3  charge  radii. 
Therefore,  to  indicate  the  adequacy  of  the  free-water  relations 
for  this  close-in  domain,  pressure-time  records  were  obtained 
from  two  experiments  conducted  in  a  closed,  rigid-wall,  cylin¬ 
drical  control  vessel  shown  in  Figure  7.8.  In  each  test: 

(a)  a  1/4"  diameter  tourmaline  piezoelectric  gage  was 
located  at  the  mld-merldlan  approximately  1/4"  from  the  internal 
wall  surface 

(b)  the  vessel  was  filled  completely  with  water  at 
atmospheric  pressure  and  70®F 

(c)  a  compact  l8-gm  pentollte  charge  was  detonated  at 
the  vessel  centroid. 
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A  pressure-time  trace  obtained  from  one  of  these  tests  Is 
shown  In  Figure  7. 9*  The  peak  shock  pressure  of  ~  52,000  psl 
was  a  combination  of  the  Incident  and  reflected  waves.  After 
the  shock  decayed,  two  steps  of  pressure  were  noted.  The  first 
step,  ~  25,000  psl.  Is  believed  to  be  the  maxlmvim  Internal- 
blast  pressure;  It  should  have  remained  at  this  level  for  the 
duration  of  the  1/2 -msec  trace*.  However,  due  to  a  partial 
top-closure  failure  and  subsequent  movement,  the  pressure 
suddenly  dropped  to  19,000  psl.  The  piezoelectric  gages 
employed  In  this  study  were  calibrated  via  the  method  and 
rationale  reported  by  R,  L,  Davis  In  reference  (s). 

Prom  a  reactor  containment  point  of  view,  the  most  Important 
portion  of  the  plot  shown  In  Figure  7.9  is  the  shock  wave. 
Conventional  free-water  equations  for  the  pentollte  shock  wave 
are  given  by  Cole  In  reference  (l)  and  by  the  Naval  Ordnance 
Laboratory  In  reference  (t)  as 


(7.1) 


where 


(7.2) 


(7.3) 


^  It  Is  to  be  expected  that  the  Internal  blast  pressure  In 
the  rigid  wall  vessel  would  exceed  significantly  that  correlative 
with  gross  plastic  deformation.  For  example.  If  the  product 
gases  In  this  rigid  wall  vessel  were  allowed  to  expand  iso- 
thermally  100  fold,  as  In  the  case  of  gross  plastic  deformation, 
the  blast  pressure  would  be  reduced  to  about  250  psl  -  several 
orders  of  magnitude  below  the  shock  pressure. 
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and  the  nomenclature  and  units  are 

P  ....  incident  wave  pressure,  psig 

Pjji  ....  peak  pressure  of  incident  wave,  psig 

t  ....  time,  msec 

9  ....  time  constant  of  incident  wave  decay,  msec 

W  ....  charge  weight,  lb 

Rj  ....  distance  from  charge,  ft 

If  the  shock  is  assumed  to  be  a  plane  acoustic  wave  reflected 
from  a  normal  rigid  surface  with  a  velocity  of  about  5OOO  ft/sec, 
then,  for  W  «  18  gm  «  0.04  lb,  the  calculated  pressure  felt  by  a 
gage  1/4”  from  a  2,5"  radius  cylinder  wall  would  be  that  given 
by  the  dashed-line  curve  shown  in  Figure  7.10. 

The  frequency  response  of  the  recording  instrumentation  was 
about  50,000  cycles/second.  Modification  of  the  calculated  curve 
was  made  to  account  for  this  recording  deficiency  (this  was  more 
feasible  than  correcting  the  experimental  data),  and  the  solid- 
line  curve  in  Figure  7.10  resulted,  A  comparison  of  this 
hypothesized  solid -line  curve  with  the  shock  portion  of  the 
actual  pressure-time  trace  (Figure  7.9)  is  significant:  good 
agreement  exists  in  the  sense  that  the  parameter  of  Importance 
to  containment  is  the  Impulse  or  the  area  under  the  incident 
pressure-time  curve.  These  experiments  provided  Important 
additional  evidence  that  conventional  free-water  shock  relations 
can  be  employed  fruitfully  in  characterizing  the  shock  loading 
function  correlative  with  containment. 

The  application  of  equations  (7Vl),  (7.2),  and  (7.3)  to 
the  vessel  and  test  conditions  specified  in  Figure  6.2  yields 
the  pressure-time  plot  shown  in  Figure  7.11.  Superposed  on 
this  curve  are  the  deformation-  and  velocity-time  plots  given 
by  Figures  6.3  and  6.4.  It  is  noted  that  the  displacement  and 
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preesure  curves  are  shown  with  the  same  origin,  althovigh  it 
is  known  that,  time -wise,  the  shock  wave  preceded  the  plastic 
deformation.  Since  the  framing  rate  of  the  hi^ -speed  camera 
could  not  define  time  zero  closer  than  40  microseconds,  common 
origins  were  chosen  for  convenience.  It  is  seen  that  90  per  cent 
of  the  incident  impulse  occurred  within  100  microseconds  as  com¬ 
pared  to  the  total  deformation  time  of  600  microseconds.  The 
fact  that  the  loading  time  was  short  relative  to  the  deformation 
time  is  consistent  with  the  instantaneous  velocity  hypothesis 
of  section  6.2:  when  the  shock  wave  strikes  the  internal  wall 
of  the  confining  vessel,  the  wall  responds  with  a  virtually 
Instauitaneous  initial  velocity. 
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FIG.  7.1  TYPICAL  TEST  VESSELS  OF  SEAMLESS  AND  WELDED  CONSTRUCTIONS 
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10"  ID  -  20"  LENGTH  -  1/4"  WALL  10"  ID  -  20"  LENGTH  -  1/4"  WALL 

304  STAINLESS  STEEL  304  STAINLESS  STEEL 

WATER  -  FI  LLED  WATER  -  FI  LLED 

CHARGE  AT  CENTROID  CHARGE  AT  CENTROID 
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10"  ID  -  20"  LENGTH  -  1/8"  WALL 
304  STAINLESS  STEEL, WELDED  PLATE 
WATER-FILLED,  CHARGE  AT  CENTROID 
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FIG.  7.3  UNIFORM  DEFORMATION  INCREASE  WITH  CHARGE  WEIGHT  FOR  WATER-FILLED  VESSELS 
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FIG.  7.4  TEST  VESSEL  FITTED  WITH  RIGID  END  CLOSURES 


NOTE:  NOZZLE  WALL  THICKNESS 
EQUAL  TO  THAT  OF  VESSEL 
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FIG,  7.5  VESSELS  FITTED  WITH  APERTURE  AND  NOZZLE 
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FIG.  7.6  EFFECTS  OF  APERTURE  AND  NOZZLE  ON  DEFORMATION 


VESSELS  SIZE:  5"  ID  -  10"  LENGTH  -  1/8"  WALL 
MATERIAL:  304  STAINLESS  STEEL,  CAST 
CHARGE  LOCATION:  CENTROID  OF  VESSEL 
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FIG.  7.7  UNIFORM  DEFORMATION  INCREASE  WITH  CHARGE  WEIGHT  FOR  AIR-FILLED  VESSELS 
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PRESSURE  (PSI  X  10 
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TIME  (fj  SEC  X  10"^) 


FIG.  7.9  RECORPED  PRESSURE-TIME  TRACE  FROM  18-Gm  EXPLOSION  IN  CONTROL  VESSEL 
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TIME  MEASURED  FROM  ARRIVAL  OF  INCIDENT  WAVE  (m  SEC) 


FIG.  7.10  PREDICTED  PRESSURE-TIME  CURVES  FOR  18-Gm  EXPLOSION  IN  CONTROL  VESSEL 
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TABLE  7.1 


RESULTS  Ol^  EXPERIMENTS  IN  EXPLOSION-LOADED,  WATER-FILLED  IDEALIZED  VESSELS 


INTERNAL 

LENGTH 

(INCHES) 

WALL 

PENTOLITE 

TERMINAL 

FAILURE 

INITIAL 

DIAMETER 

THICKNESS 

CHARGE 

RADIAL 

OR 

STRAIN 

(INCHES) 

(INCHES) 

WEIGHT 

STRAIN 

CONTAINMENT 

RATE 

(GRAMS) 

(IN/IN) 

MODE 

(SEC-^) 

304  STAINLESS  STEEL,  WELDED-PLATE  CYLINDERS 
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0.289 

5 

- 

1 

0.316 

1/3 

2070 

2 

0.343 

1,4 

- 

3 

0.366 

1,4 

- 

4 

0.195 

6 

- 

5 

0.271 

6 

- 

6 

0.331 

6 

- 

7 

0.349 

6 

- 

8 

0.357 

5 

2850 

9 

0.304 

1,3 

- 

10 

0.349 

1,4 

- 

11 

0.295 

6 

. 

12 

0.338 

5 

2900 

13 

0.295 

1,4 

- 

14 

0.338 

1,4 

- 

15 

0.366 

6 

- 

16 

0.434 

6 

- 

17 

0.506 

6 

1570 

18 

0.568 

5 

- 

19 

0.280 

6 

. 

20 

0.315 

6 

- 

21 

0.348 

6 

1300 

22 

0.383 

6 

- 

23 

0.425 

6 

- 

24 

0.487 

5 

1570 

25 

0.387 

2,4 

1870 

26 

0.546 

2,4 

- 

27 

0.569 

2,4 

1940 

28 

0.512 

5 

1640 

29 

0.468 

3 

- 

30 

0.490 

4 

1970 

31 

0.180 

6 

32 

0.231 

2,3 

568 

33 

- 

2,4 

790 

34 

0.297 

2,4 

805 

35 

0.412 

6 

_ 

36 

0.523 

5 

- 

37 

0.516 

3 

38 

0.529 

5 

39 

0.561 

4 

- 

40 

0.110 

6 

41 

0.265 

5 

- 

42 

0.288 

2,3 

- 

43 

0.342 

2,4 

- 

44 

0.236 

1,3 

- 

45 
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TABLE  7.1  (Continued) 


RESULTS  OF  EXPERIMENTS  IN  EX  PLOSION- LOADED,  WATER-FILLED  IDEALIZED  VESSELS 


INTERNAL 
DIAMETER 
( INCHES) 


WALL 

PENTOLITE 

TERMINAL 

FAILURE 

INITIAL 

LENGTH 

THICKNESS 

CHARGE 

RADIAL 

OR 

STRAIN 

(INCHES) 

(INCHES) 

WEIGHT 

STRAIN 

CONTAINMENT 

RATE  , 

(GRAMS) 

(IN/IN) 

MODE 

(SEC*') 

304  STAINLESS  STEEL,  CENTRIFUGALLY-CAST  CYLINDERS 


TEST 


NO. 


10 


- 0.125 


- 0.250 


20 - 0.125 


-  0.250 


10 - 20 


-  0.500 


16.0 

0.295 

6 

- 

46 

18.0 

0.352 

5 

- 

47 

19.0 

0.381 

3 

- 

48 

21.0 

- 

4 

- 

49 

50.1 

0.375 

5 

_ 

50 

50.5 

0.354 

5 

- 

51 

50.5  *** 

0.385 

3 

- 

52 

58.1 

0.425 

6 

- 

53 

70.2 

0.461 

6 

- 

54 

75.3 

0.497 

5 

- 

55 

85.0 

0.483 

3 

- 

56 

_  100.0 

- 

4 

- 

57 

Q  19.0 

0.370 

5 

- 

58 

’*189.8 

0.463 

5 

1450 

59 

189.9 

0.440 

3 

1490 

60 

200.2 

0.413 

4 

- 

61 

—250.1 

- 

4 

- 

62 

“461.0 

0.421 

6 

- 

63 

481.2 

0.447 

5 

1810 

64 

497.7 

0.436 

3 

- 

65 

569.2 

0.440 

4 

- 

66 

410  STAINLESS  STEEL,  CENTRIFUGALLY-CAST  CYLINDERS 
ANNEALED 


5 


10 


10 - 0.250 


20 


1—0.250 


*—0.500 


25.3 

0.118 

6 

- 

67 

28.0 

0.131 

6 

- 

68 

28.0 

0.132 

6 

1860 

69 

28.0 

0.141 

5 

- 

70 

29.0 

- 

3 

- 

71 

,  30.3 

- 

4 

- 

72 

110.2 

0.173 

6 

- 

73 

130.1 

0.220 

6 

- 

74 

160.6 

0.258 

6 

310 

75 

180.4 

0.292 

6 

- 

76 

199.6 

0.341 

6 

- 

77 

220.5 

0.364 

5 

- 

78 

239.8 

0.379 

3 

- 

79 

‘210.0 

0.125 

6 

- 

80 

230.4 

0.147 

5 

950 

81 

249.7 

. 

4 

- 

82 

320.5 

- 

4 

- 

83 

HARDENED 


5 - 10 


0.250 


22.0 

26.0 


0.045 

0.050 


5 

3 


84 

85 
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TABLE  7.1  (Continued) 

RESULTS  OF  EXPERIMENTS  IN  EXPLOSION-LOADED,  WATER-FILLED  IDEALIZED  VESSELS 


INTERNAL 

DIAMETER 

(INCHES) 


10 


20 


10 


15 


NOTES 


LENGTH 
( INCHES) 


WALL 

THICKNESS 

(INCHES) 


20 


0.515  — 


1—0.515  — 


40- 


L-  1.030 


MILD  STEEL,  WELDED  PLATE  CYLINDERS 

0.271 


20 


0.253- 


140.0  ♦♦ 
145.0  ♦♦ 
145.0 


1.3 

1.4 
1,4 


6061  -  T6  ALUMINUM,  EXTRUDED  CYLINDERS 


30 


0.500- 


60.0  ** 
85.2 
90.0  ** 
99.9 


0.0553 

0.0855 

0.0875 

0.0929 


6 

5 

3 

4 


1010 


*  Cylinders  were  of  double-weld  construction  (two  longitudinol  seoms,  180  opart). 
**  No  certified  mechonicol  properties  were  avoiloble  for  this  materiol, 

*'**  Cylinders  were  fitted  with  rigid,  radiol  end  closures. 

+  The  terminol  strain  is  the  ot-rest  mid-meridian  strain. 


PENTOLITE 

TERMINAL 

FAILURE 

INITIAL 

CHARGE 

RADIAL 

OR 

STRAIN 

TEST 

WEIGHT 

STRAIN 

CONTAINMENT 

RATE  , 
(SEC“') 

NO. 

(GRAMS) 

(IN/IN) 

MODE 

WELDED  PLATE  CYLINDERS 

200.0  ♦ 

0.154 

6 

- 

86 

229.3  ♦ 

0.183 

5 

867 

87 

275.0  * 

0.118 

1,3 

- 

88 

__324.4  ♦ 

0.128 

1,4 

- 

89 

■“902.7 

0.242 

6 

- 

90 

1103.0 

0.300 

6 

- 

91 

1200.3 

0.328 

5 

670 

92 

J351.5 

0.357 

3 

650 

93 

"2002.6 

- 

1,4 

- 

94 

2251.0 

0.230 

5 

- 

95 

2497.0 

0.226 

1,3 

- 

96 

2498.0 

0.185 

2,3 

- 

97 

98 

99 

100 


101 

102 

103 

104 


++  The  following  code  describes  mode  of  cylinder  failure  or  contoinment. 

1  - failure  occurred  in  heat  effected  zone  adjocent  to  weld 

2  ---  foilure  occurred  in  weld  material 

3  ---  morginol  rupture 

4  ---  non-morginol  rupture 

5  — -  marginol  contoinment 

6  ---  non-marginol  contoinment 
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.  All  charges  were  pressed  pentalite. 
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Vessel  material  was  304  stainless  steel,  cast. 

These  experiments  were  conducted  In  vessels  of  the 
type  shown  In  Figure  7.4. 

All  charges  were  pressed  pentolite. 
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CHAPTER  8 


EVALUATION  OP  EMBRYONIC  FUNCTIONS 


The  extensive  experimental  results  of  Chapter  7  provide 
the  necessary  and  sufficient  botmdary  conditions  for  evaluating 
the  erabiTTonlc  f\mctlona,  and  't'2^  generated  via  Salnt-Venant  ’ s 
semi-inverse  method  in  Chapter  6. 

8,1  Efficiency  Factor  Function  ♦g.  It  proves  convenient  to 
evaluate  the  embryonic  functions  in  the  reverse  order  of  their 
derivation.  We  start  by  rewriting  equation  (6.l8)  in  the  fonn 


„4/3 


(8.1) 


*1*2  “  «’t  '("i  -  •’i)'”' 


where  the  true  stress  function  as  previously  defined,  is 


(6.14) 


All  parametric  values  on  the  right-hand  side  of  equation  (8,1) 
are  known  from  Table  7.1  for  the  case  of  containment,  and  the 
calculated  products  these  cases  are  given  in  Table  C-1 

of  Appendix  C. 

It  was  hypothesized  in  equation  (6.19)  that  ^2  ^  function 

of  Ri/lio*  flg'ire  8.1  shows  this  significant  dependence,  from 
which  the  best  straight  line  through  averaged  values  of  the 
groupings  is  taken  to  be 


♦2  -  10"^  (3.41  +  0.117  RiAq) 


(8.2) 
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8.2  Deformation  Energy  Function  ♦ i .  For  cylinders,  we  define 
the  radial  strain  e  as 


.  ,  R  -  Rl 


(8.3) 


Where  R  and  are  the  Instantaneous  and  Initial  Internal 
radii,  respectively.  Then  the  strain  rate  t  is 


1  1  dR  V 

®  “  RT  “  iq’ 

where  V  Is  the  scalar  quantity,  v  «  |v|  .  Since  the  Initial 
strain  rate  of  the  wall  at  the  mld-merldlan  was  hypothesized 
to  Index  the  density  of  the  deformation  energy ,  we  put 


®o  "  Vo(0)/Ri  ,  (8.^) 

Combining  equations  (6.7)>  (6.8),  (6.13),  and  (8.4), 
taking  K  «  2.18  for  pentollte  from  reference  (l),  and  account¬ 
ing  for  units,  we  obtain 


,  (8.5) 

Rl(Re  -  Rl) 
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Table  C-2  In  Appendix  C  gives  5  for  23  experiments  In  which 
the  strain  rate  was  monitored.  The  trial  format  assumed 
?  =  ?  (Rl/ho)«  This  functional  relation  Is  plotted  In  Figure 
8.2;  the  best  straight  line  through  the  averaged  values  of  the 
groupings  yields 


?  -  1.47  +  0.0373  RiAo 


(8.6) 


and  upon  combination  with  (8.5),  we  find 


_ 20216 

(1.47  +  0.0373  RiAo)  w  %(Re  ■  ^1  ^ 


(8.7) 


The  strain  rate  based  on  (8.7)  appears  In  Table  C-1  for  all 
contained  experiments.  Table  C-2  gives  the  deviation  of  these 
(8.7)  results  from  the  experimental  values,  and  It  Is  noted 
that  the  deviations  are  as  high  as  15  per  cent,  A  plot  of 
and  both  from  Table  C-1,  Is  given  In  Figure  8.3;  motivation 
for  the  log-log  scales  was  provided  by  the  fractional-power 
strain  rate  function  given  In  equation  (6.16).  The  best  straight 
line  throvigh  these  data  yields 


-  0.318 


(8.8) 


Which,  when  Incorporated  Into  (8.I),  reduces  the  effect  of  the 
maxlmxim  Strain  rate  deviation  of  (8.7)  from  15  per  cent  to  less 
than  3  per  cent. 
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The  simultaneous  solution  of  (8.7)  and  (8.8)  yields 


♦l 


0.318 


-1  - 

ITH^). 


0.15 


(8.9) 


and  we  now  have  and  ^2  completely  defined  In  terms  of  charge 
weight  and  simple  properties  of  the  confining  vessel. 
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FIG.  8.1  VARIATION  OF  W,  %  WITH  R./h 
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FIG.  8.2  VARIATION  OF  ^  WITH  R./h 
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FIG.  8.3  VARIATION  OF  WITH  INITIAL  STRAIN  RATE  £ 
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CHAPTER  9 
CONTAINMENT  LAWS 

Th»  TNT  rationale,  embryonic  format,  extensive  experi¬ 
mental  program,  and  evaluation  of  the  trial  functions  now 
culminate  In  final  Containment  Law  formulations. 

9.1  TNT  Containment  Law  for  Ideal  Vessels.  The  embryonic 
form  of  the  Containment  Law 


W  - 


(6.18) 


constitutes  certain  necessary  conditions;  further  necessary 
conditions  are  provided  by 


.?7^+  ^0®0'373^ii/ho-)  ] 


0.15 


(8.9) 


♦2  •=  10"^  (3.41  +  0.117  RiAo) 


(8.2) 


The  simultaneous  solution  of  equations  (6.I8),  (8.2),  and  (8.9) 
yields  the  TNT  Containment  Law  for  ideal  vessels  (no  nozzles, 
weldments,  or  rigid  end  constraints): 


1.407  Of:  «  (3.41  4  0.117  Ri/ho)(Rl  - 

10^  (1.47  +  0.0373  Ri/ho)^^^Ri^*^  - 


0.811 


(9.1) 
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where 

«t  -  "y  +  [■’u  (1  -  oy 


and  the  bar  on  W  Indexes  the  valid  domain  of  the  equation  as 
set  forth  In  (a)  through  (g)  of  section  6,3.  Theoretically, 
the  maxlm\im  value  of  (9.1)  occurs  for  ®  *  ®u  given  by 


l.»0T  Oa  (1  0.117  V^oXl^  - 

ID^  (1.^7  ♦  0.0373 


0.611 


(9.2) 


In  practice,  however,  the  mid -meridian  strain  correlative  with 
marginal  containment  Is  usually  less  than  ®u.  This  Is  seen 
from  comparing  the  conventional  ultimate  strain  values  of 
section  7.1  with  the  marginal  containment  strains  reported  In 
Table  7.1 

9.2  TNT  Containment  Law  for  Real  Vessels.  For  a  containment 
law  to  be  useful  to  the  pressure  vessel  designer.  It  must  yield 
safe  solutions  for  "real"  vessels,  l.e.,  vessels  fitted  with 
nozzles,  weldments,  and  rigid  end  closures.  Experiments 
designed  to  evaluate  specifically  the  limiting  effects  of 
these  discontinuities  on  containment  were  reported  In  Chapter  7. 
On  the  basis  of  the  results.  It  was  found  that  the  safe  or 
permissible  strain  domains  correlative  with  the  discontinuities 
are 
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(a)  waldftd  nozzles:  «  <  1/2  (section  7.4) 

(b)  rigid  end  closures:  1/3  ®u  ^  ®  ^  1/2  *u  (Bectlon  7.3) 

(c)  welded  plate:  «  <  1/3  Su  (zactlon  7.2) 

From  an  examination  of  these  limits.  It  Is  postulated  that 
real  vessels  can  be  expected  to  withstand  safely  strains  as 
large  as  1/3  ®u.  ^posing  this  restriction  on  equation  (9*1), 
we  find  that  the  TMT  Containment  Law  for  real  vessels  Is 

-  Wj  (9.3) 


except  that 

(a)  e  <  1/3  _ 

(b)  the  valid  domain  Indexed  by  the  bar  on  consists 
of  (a)  through  (f)  of  section  6.3,  and  additionally 

(c)  the  vessel  wall,  weldments,  closures,  and  nozzles 
must  possess  everywhere  the  chemical,  mechanical,  and  physical 
properties  specified  In  the  design,  and 

(d)  these  properties,  as  they  relate  to  containment,  must 
not  be  negated  In  time  by  flaws,  NDT  (nil  ductility  temperature) 
considerations,  radiation  effects,  stress  concentrations,  high 
temperatures,  cyclliig,  etc. 

The  literature  abovinds  with  studies  on  phenomena  that 
degrade  the  mechanical  properties  of  steel.  Among  the  more 
authoritative  of  these  are  0,  R.  Irwin's  treatment  of  flaws  In 
reference  (u),  the  work  of  W,  S,  Pelllnl  and  P.  P.  Puzak  on 
NDT  considerations  In  reference  (v) ,  and  Q,  M.  Adamson's  data 
on  radiation  effects  In  reference  (w). 
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The  maximum  value  of  equation  (9.3)  occurs  for  e  =  1/3  ®u 
and  Is  given  by 


0.1563  ♦  0.117 

lo5  (2  oy  +  0.0373 


0.611 

(9.4) 


where.  In  review,  the  nomenclature  and  units  for  equations  (9.l) 
through  (9.4)  are 

W  ....  charge  weight  (TNT  or  pentollte),  lb 

w  ....  weight  density  of  vessel  material,  Ib/ft^ 

Ri  ....  Initial  Internal  radius  of  vessel,  ft 

Rg  ....  Initial  external  radius  of  vessel,  ft 

ho  ....  Initial  wall  thickness  of  vessel,  ft 

c  ....  conventional*  permissible  strain  of  vessel 

material,  In/ln 

®u  ....  conventional*  ultimate  strain  of  vessel  material, 
In/ln 

or^  ....  true  stress*,  psl,  see  equation  (9.1) 

Cy  ....  conventional*  yield  stress  of  vessel  material,  psl 

....  conventional*  ultimate  stress  of  vessel  material, 
psl 

*  It  Is  noted  that  all  Ingredients  of  tlie  Containment  Laws  are 
In  terms  of  conventional  materials  and  properties  -  quantities 
readily  available  to  the  design  engineer;  however.  If  the  vessel 
material  Is  subjected  to  temperatures  that  degraule  Its  energy 
absorption  properties  -  l.e..  Its  stress-strain  values  -  then, 
of  course,  these  degraded  values  must  be  employed  In  equations 
(9.1)  throxjgh  (9.4). 
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9.3  Accuracy.  Deviations  between  predictions  based  on  equation 

(9.1)  and  results  from  the  contained  experiments  are  given  In 

Table  C-1  aixi  plotted  In  Figure  9*1*  If  all  tests  are  considered, 

the  average  deviation  Is  ^  per  cent,  and  the  maximum  deviation 
-♦■14 

Is  _2o  cent.  If  we  exclude  those  tests  for  which  no  certified 
mechanical  properties  were  available  (denoted  by  circles)  and 
those  tests  for  which  the  final  at-rest  strain  exceeded  the 
maxlmvim  certified  ultimate  elongation  (denoted  by  sqviares),  the 
deviation  Is  substantially  reduced  to  an  average  +  4  per  cent 
and  a  maxvimum  of  P®^  cent.  Irrespective  of  the  preferred 
deviation,  agreement  between  the  containment  law  predictions 
and  the  experimental  results  Is  considered  excellent  In  light 
of  the  following  statements : 

(a)  The  calculations  made  In  Table  C-1  of  Appendix  C  were 
based  on  an  average  stress-strain  equation  for  a  given  material. 
Frequently,  Individual  certified  tensile  tests  varied  from  the 
average  by  as  much  as  10  per  cent  and  In  a  few  cases  as  high  as 
20  per  cent. 

(b)  Deviations  In  the  replicability  of  explosion  phenomena 
frequently  exceed  5  per  cent  and  are  sometimes  as  high  as  10  per 
cent. 

Experimental  results  of  vessels  fitted  with  weldments  and 
rigid  end  closures  are  Included  in  Table  C-1  and  Figure  9*1  and, 
hence.  In  the  accuracy  data  presented  above.  Experimental 
results  of  the  vessels  fitted  with  nosszles  also  satisfy  equation 

(9.1)  with  a  maximum  deviation  of  1  per  cent.  A  typical  example 
Is  taken  from  Table  7.2:  304  stainless  steel  (cast)  vessel, 

10"  ID,  20"  length,  1/4"  wall,  e  «  0.328  In/ln,  W  -  I30  gm. 

With  o^  determined  from  values  given  In  Appendix  C,  equation 

(9.1)  yields  Wj  «  129  gm,  a  deviation  of  less  than  1  per  cent. 

The  strain  profile  of  this  vessel  Is  shown  In  Figure  7.6;  since 
lx;  Is  asymetrlc,  c  indexes  the  "local"  or  uninhibited  portion 
of  the  mld-merldlan  (which.  In  this  case,  lies  opposite  the  one 
nozzle).  With  two  or  more  nozzles,  the  local  region  would  be 
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more  restricted.  Because  of  this  difference  in  the  interpreta¬ 
tion  of  c,  the  deviation  data  for  nozzle  experiments  were  not 
included  in  Table  C-1  and  Figure  9*1 

9.4  Upper-Bound  Containment  Law  for  Reactor  Vessels.  It  is 
recalled  that  the  excursion  accident  spectmam  is  large  in  both 
magnitude  and  flxix  of  the  energy  released;  it  ranges  from 
small  accidents  to  thousands  of  megawatt-seconds  occurring  in 
times  from  a  few  microseconds  to  seconds  or  more,  and  these  can 
take  place  anywhere  in  the  coolant-loss  domain  from  no  loss  to 
a  complete  loss.  The  range  of  excursions  and  the  conditions 
under  which  they  may  occur,  the  multiplicity  and  complexity  of 
the  parameters,  and  the  state  of  the  art  preclude  an  exact 
solution  to  the  reactor  containment  problem  at  this  time. 

This  means  that  only  restricted  solutions  are  possible. 

Of  these,  we  can  postulate,  a  priori,  that  upper-bound  solutions 
for  containing  the  postulated  excursion  under  the  most  adverse 
conditions : 

(a)  energy  flux  comparable  to  that  of  TNT 

(b)  reactor  vessel  filled  with  liquid  coolant 

are  the  most  basic  and  the  most  valuable;  and  that  these 
solutions  are  precisely  those  provided  by  the  TNT  Containment 
Laws,  (9.1)  through  (9.4).  One  further  clarification  is 
necessary; 

To  avoid  ambiguity  in  relating  a  postulated 
excursion  to  a  TNT  simulation,  it  is  assumed 
that,  irrespective  of  the  compounding  of 
nuclear  or  other  events  that  may  occur  in 
the  reactor,  the  upper  bound  on  the  resultant 
total  accident  (the  energy  released)  is  in 
each  case  relative  to  and  occurs  at  the  normal 
operating  conditions  of  the  reactor.  Further¬ 
more,  it  is  assumed  that  the  release  of  gas 
products  per  unit  energy  of  reactant  does  not 
exceed  that  for  TNT. 
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Upon  some  reflection.  It  is  seen  that  the  limitations  on 
compovindlng  and  gas  products  are  necessary  for  upper-bound 
containment  predictions;  but  that  since  the  majority  of 
excursions  can  be  Interpreted  to  satisfy  these  conditions, 
no  serious  loss  In  generality  accrues.  It  Is  concluded  that 
the  TNT  Containment  Laws  (9.3)  and  (9.4)  constitute  a  safe 
and  wholly  defensible  upper-bound  solution  to  the  entire 
spectrum  of  postulated  excursions  and  the  conditions  xinder 
which  they  may  occur. 

9.5  Less  Adverse  Accident  Conditions.  The  TNT  Containment 
Laws  (9.3)  and  (9.4)  will  provide  safe  upper  bounds  for  the 
entire  spectmm  of  excursions  and  the  conditions  under  which 
they  may  occur.  If  the  energy  flux  Is  similar  to  that  of  TNT 
and  the  reactor  vessel  Is  full  of  liquid  coolant,  these  Laws 
will  yield  upper  bovinds  that  are  both  safe  and  reasonable. 

If,  however,  the  excursion  conditions  are  less  adverse  -  the 
energy  Is  released  at  rates  less  than  that  of  TNT  and/or  a 
loss  of  coolant  occurs  -  the  upper-bound  solution  will,  of 
course,  be  safe,  but  Its  reasonableness  will  be  decreased. 
Modifications  of  the  upper-bound  Laws  are  necessary  If  solu¬ 
tions  for  less  adverse  conditions  are  to  be  both  safe  and 
reasonable.  Three  cases  of  excursion  conditions  less  adverse 
than  those  for  the  upper-bound  solution  are  presented  In  the 
following  sections. 

9.6  TNT  Energy  Flux,  Complete  Loss  of  Coolant.  The  explosion 
experiments  conducted  In  vessels  filled  with  air  Instead  of 
water  (section  7.5)  provide  pertinent  Information  for  the  case 
of  TNT  energy  flux  and  a  complete  loss  of  coolant  from  the 
reactor  vessel  prior  to  the  accident.  A  direct  comparison 
between  the  water  and  air  conditions  Is  presented  In  Figure  9.2. 
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It  is  assumed  here  that  the  pressure  within  the  reactor  vessel 
after  a  complete  loss  of  coolant,  but  prior  to  an  excursion, 
is  atmospheric.  This  assumption  is  not  restrictive  since  only 
sodium  reactors  are  subject  to  postulated  excursions  following 
a  complete  loss  of  coolant,  in  which  case  the  pressure  would  be 
at  or  near  atmospheric.  It  is  seen  Immediately  that  everywhere 
within  the  specified  radial  strain  domain  for  air,  the  contain¬ 
ment  potential  is  at  least  a  factor  of  three  greater  than  that 
for  water  and  that  up  to  0.20  in/in,  the  factor  is  at  least 
four.  These  values  for  ideal  vessels,  however,  could  not  be 
expected  to  hold  for  real  vessels. 

Since  the  localization  factor  for  air  is  half  that  of 
water  (Figure  6.1),  deformation  at  the  mid -meridian  would  be 
much  more  violent;  hence,  nozzle  and  weld  failures  would  occur 
at  smaller  radial  strains  for  air  than  for  water.  This  effect 
can  be  compensated  for  by  restricting  the  strain  for  the  air 
case  to  e  <  ®\i/6.  Then  the  permissible  values  for  304  stainless 
steel  are 

®alr  =  ^  ®water  =  ®u/3  =  0.20  in/ln 

and  from  a  comparison  of  the  correlative  charge  weights  (Figure 
9.2),  it  is  hypothesized  that  the  value  2  constitutes  a  safe 
upper  bound  everywhere  in  the  radial  strain  domain  of  equation 
(9.3) ^  and,  of  course,  this  includes  (9.4). 

The  alr/water  containment  ratio  decreases  with  increasing 
radial  strain  as  shown  in  Figure  9.2,  and  beyond  the  limiting 
strain  for  air  (about  0.2?  in/ln),  the  decrease  becomes  signifi¬ 
cant.  For  example,  a  typical  limiting  strain  for  water  (in 
the  case  of  no  rigid  end  constraints)  is  0.46  in/in,  from  which 
the  Figure  yields  an  air/water  containment  ratio  of  about  1.9. 
Such  values  are  somewhat  academic,  however,  since  they  are  not 
consistent  with  admissible  strains  of  the  realistic  Containment 
Law,  (9.3). 
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For  Intermediate  losses  of  liquid  coolant,  the  composite- 
vessel  experiments  of  reference  (1)  Indicate  that  the  condition 
of  overriding  Importance  Is  whether  or  not  the  coolant  level 
In  the  vdssel  Is  sufficient  to  completely  envelope  the  core. 

If  the  core  Is  completely  enveloped,  the  containment  potential 
will  follow  closely  that  of  the  no  coolant-loss  condition  with, 
perhaps,  a  slight  enhancement.  If  the  coolant  loss  Is  such 
that  the  core  Is  not  enveloped,  the  containment  potential  will 
be  somewhat  greater  (say  3/2)  than  that  for  the  no-loss  condi¬ 
tion,  and  this  potential  factor  will  converge  to  the  value  2 
for  Increasing  coolant  loss. 

In  the  containment  domain,  the  strength  of  the  shock  Is  not 
an  appreciable  function  of  the  temperature  of  the  medium  through 
which  It  propagates.  In  view  of  this.  It  Is  postulated  (for 
TNT  energy  flux)  that  the  values  for  no  coolant  loss  and  2 
for  complete  coolant  loss  are  valid  everywhere  In  the  temperature 
range  from  65®P  to  850®P. 

9.7  Lesser  Energy  Flux,  No  Coolant  Loss.  Experiments  with 
solid  propellant  In  water-filled,  5”  ID  vessels  provide  useful 
bounds  for  the  case  of  energy  fluxes  less  than  that  of  TNT  and 
no  loss  of  liquid  coolant  from  the  vessel  prior  to  the  accident. 
The  propellant  employed  In  the  experiments  of  reference  (k) 
possessed  an  energy  density  similar  to  that  of  TNT;  the  water 
was  70*^.  The  energy  flux  Is  reflected  In  the  pressure-time 
and  strain- time  plots  of  Figure  9.3.  The  release  time  of  about 
40  milliseconds.  If  scaled  to  a  typical  reactor  -  say  one  with 
ID  of  10  feet  -  would  correspond  to  about  one  second.  More 
generally,  for  the  conditions  here,  the  release  time  for  a 
reactor  of  radius  Rj^  (feet)  can  be  taken  as 


«=  200  (tnllllseconds ) 
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where  la  the  numerical  value  of  Rj .  The  propellant 
experiments  and  the  TNT  Containment  Law  for  Ideal  Vessels, 
equation  (9.1) >  Index  In  Figure  9.4  the  significant  Increase 
In  containment  correlative  with  energy  fluxes  less  than  “the 
TNT  flux. 

Although  the  propellant  experiments  were  conducted  In 
vessels  with  rigid  end  closures,  the  uniformity  of  loading  over 
the  full  vessel  length  resulted  In  deformations  similar  to 
those  of  the  open-end  explosion  experiments.  It  Is  seen  from 
Figure  9.4  that  everywhere  within  the  specified  radial  strain 
dcmaln  for  the  propellant,  the  containment  potential  Is  at 
least  a  factor  of  five  greater  than  that  for  the  explosive. 

Much  more  Important,  however.  Is  the  restricted  strain  domain, 

•  ^  ®u/3»  TNT  Containment  Law  for  real  vessels,  equation 

(9.3).  For  304  stainless  steel,  an  allowable  maximum  Is 
«  -  0.20  In/ln.  From  the  use  of  this  value  In  comparing  the 
correlative  weights  of  propellant  and  explosive  In  Figure  9.4, 

It  Is  hypothesized,  for  water  at  70®F  and  energy  release  times 
>  that  the  value  8  1  I  Is  a  safe  upper  bound,  |  %  |  being 

the  energy  density  of  without  regard  to  energy  flxix. 

Unlike  the  case  for  TNT  energy  flux,  containment  potential 
for  lesser  fliixes  decreases  significantly  with  Increased  system 
temperature.  Liquid  coolant  temperatures  In  excess  of  70®F 
will  decrease  the  enhancement  factor  of  eight  In  accord  with 
the  basic  heat-sink  rationale  presented  In  Appendix  D  and 
reference  (l).  In  essence  this  rationale  hypothesizes  that, 
for  extensive  metal  and/or  coolant  heat  sinks,  the  Internal 
blast  pressure  varies  directly  with  absolute  temperature  and 
charge  weight.  Thus  for  water  at  200®F  and  release  times 
a  Cj,  It  Is  hypothesized  that  6  |  ^  |  Is  a  safe  upper  boxind. 

For  sodium  at  850®F  and  release  times  at  fij.  It  Is  hypothesized 
that  a  safe  upper  bound  Is  3  I  I  . 
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It  Is  recalled  that  the  principal  reason  for  considering 
the  less  adverse  conditions  of  the  accident  spectrum  was  to 
generate  containment  law  solutions  for  this  portion  of  the 
spectrum  that  would  be  more  reasonable  than  that  given  by  the 
TNT  law.  In  this  vein.  It  Is  desirable  to  entertain  briefly 
the  relation 

(  ^1  )c  -  Rln  (milliseconds) 

Now  (  )q  1*  qualitatively  more  conservative  than  Oj,  but 
It  Is  believed  that  the  use  of  these  values  In  the  foregoing 
containment  predictions  would  produce  solutions  of  comparable 
reasonableness.  This  deduction  expands  the  release  time  spectrum 
of  propellant  energy  fluxes  by  5/3  and  permits  the  convenient 
tabular  format  employed  In  Table  9.1» 

The  heat-sink  experiments  of  reference  (l)  provide  a 
significant  Index  to  the  time-heat  absorption  character  of 
accidents  that  occur  In  times  of  the  order  of  hundreds  of 
milliseconds.  The  results  of  these  experiments  In  3.5”  ID 
models  Indicate  that  the  major  portion  of  heat  absorption  from 
the  Internal  blast  pressure  gases  occurred  within  five  milli¬ 
seconds.  Scaled  to  a  full-size  reactor,  say  one  with  ID  of 
10  feet,  this  time  would  be  of  the  order  of  170  milliseconds. 

More  generally,  the  release  time  here  for  a  reactor  of  radius 
Rl  (feet)  can  be  taken  as 


02  *  Rin  (milliseconds) 
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It  appears,  then,  that  no  appreciable  quantity  of  heat  la 
absorbed  In  the  Interim  between  Oi  and  and  It  Is  hypothe¬ 
sized  that  6  1  Wp  I  Is  also  a  safe  upper  bound  for  water  at 
200®P  and  release  times  ^  02*  Likewise,  It  seems  reasonable 
to  take  3  I  Wp  I  as  a  safe  upper  bound  for  sodium  at  850®P 
and  release  times  k  02 • 

At  some  point  In  the  energy  release  time  domain  less  than 
02  ,  shock  waves  will  be  produced  rather  than  a  uniform  pressure 
bu^ld-up.  This  will  be  true  Irrespective  of  whether  the  reactor 
liquid  coolant  be  lost  or  retained  In  any  quantity.  Prom  a 
consideration  of  the  dearth  of  definitive  Information  In  this 
domain.  It  Is  postulated  that  *  Wj^  I  for  the  case  of  no  coolant 
loss  and  2  1  Wj^  |  for  the  case  of  complete  coolant  loss  are  the 
highest  safe  upper  bounds  that  can  be  adequately  defended. 

9.8  Lesser  Energy  Flux,  Complete  Loss  of  Coolant.  The  explosion 
experiments  conducted  In  vessels  filled  with  air  Instead  of  water 
(section  7.5)  also  provide  some  guide  lines  for  the  case  of 
energy  fluxes  less  than  that  of  TNT  and  a  complete  loss  of 
coolant  from  the  reactor  vessel  prior  to  the  accident.  In  this 
regard  It  Is  pertinent  to  note  that ; 

(a)  In  the  containment  domain,  the  magnitude  of  the 
Internal  blast  pressure  Is  never  sufficient  to  produce  vessel 
deformation  beyond  that  generated  by  the  shock  wave;  this  Is 
true  Irrespective  of  whether  the  transmission  medium  be  liquid 
or  gas. 

(b)  Unlike  the  case  for  TNT  energy  flux  (section  S.6), 
the  lesser  energy  flux  In  air  would  not  produce  a  highly  local¬ 
ized  deformation  at  the  mid -meridian,  but  rather  a  pattern 
similar  to  that  for  TNT  In  water  (Figure  6.I).  Hence  the 
permissible  strain  for  the  lesser  energy  flxxx  In  air  and  the 
TNT  flux  In  water  would  be  the  same:  e  ^  ®u/3. 
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In  view  of  (a)  and  (b).  It  is  hypothesized  froin  Figure  9.2 
that  3  I  Wp  I  Is  a  safe  upper  bound  for  air  at  70®P  and  release 
times  i  02*  Again  as  In  section  9.6,  It  Is  assumed  that 
atmospheric  pressure  exists  In  the  reactor  vessel  prior  to 
the  excursion. 

If  extensive  metallic  heat  sinks  are  available  to  the 
released  gases,  as  In  the  case  with  the  core  region  of  the 
Enrico  Fermi  Reactor  (compare  Tables  8.2  and  8.5  of  reference 
(l)),  and  If  It  Is  assumed  that  the  partition  of  energy  with 
respect  to  the  heat  sink  for  the  slow  release  Is  similar  to 
that  for  TNT,  the  containment  potential  could  be  Increased  by 
a  factor  exceeding  two.  Under  these  conditions  It  Is  hypothe¬ 
sized  that  6  I  Wp  I  Is  a  safe  upper  bound  for  air  at  70®F  and 
release  times  i  ^2*  As  mentioned  previously,  containment 
potential  Is  affected  adversely  for  the  case  of  slow  energy 
fluxes  (those  less  than  that  of  TNT)  occurring  at  elevated 
temperatures.  The  decrease  In  containment  potential  with 
Increasing  temperature  Is  significant  Irrespective  of  whether 
the  coolant  be  retained  or  lost.  For  air  at  200®F,  extensive 
metallic  heat  sinks,  and  release  times  i  hypothesized 

again  from  the  heat-sink  rationale  of  Appendix  D  -  that  4  |  | 

Is  a  safe  upper  bound.  Similarly,  3  I  Wr  |  Is  hypothesized  to 
be  a  safe  upper  bound  for  air  at  850®P. 

For  Intermediate  losses  of  liquid  coolant,  the  containment 
potential  would  be  analogous  to  that  for  the  TNT  energy  flux: 

(a)  If  the  quantity  of  coolant  In  the  vessel  Is  sufficient 
to  completely  submerge  the  core  and  the  release  times  are  >  (12 » 
the  containment  potential  will  converge  on  the  values  for  no 
coolant  loss  -  8  1  Wp  1  for  water  at  70®F,  6  I  Wp  1  for  water  at 
200®F,  and  3  I  Wr  )  for  sodium  at  850®F  -  as  the  coolant  loss 
decreases . 
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(b)  If  the  quantity  of  coolant  In  the  vessel  Is  not 
sufficient  to  submerge  the  core  and  the  release  times  are 
i  ^2,  the  containment  potential  will  converge  on  the  values 
for  complete  coolant  loss  -  6  |  |  for  water  at  70®P,  4  |  Wp  I 

for  water  at  200®F,  and  3  I  Wp  |  for  sodium  at  850®P  -  as  the 
coolant  loss  Increases. 

9.9  Restrictions  on  Less  Adverse  Conditions .  A  summary 
of  maximum  containable  energy  releases  for  various  accident 
conditions  Is  presented  In  Table  9.1.  Por  the  case  of  no 
coolant  loss  and  energy  release  times  n  (milliseconds)  given  by 


TNT  release  time  < 


n  < 


100 

3 


^In 


the  results  In  Table  9.I  are  both  reasonable  and  highly  defensi¬ 
ble;  they  have  been  verified  directly  by  established  rationale 
and  extensive  experimental  results.  For  all  other  conditions 

100  ^  3  n  /  Rin  ^ 


the  results  In  Table  9.1  are  also  believed  to  be  safe  and 
reasonable.  However,  they  are  not  postulated  results;  they 
are  hypothesized  results.  As  such,  they  do  not  constitute 
absolute  upper  bounds  In  the  sense  of  the  postulated  quanti¬ 
ties  nor  do  they  possess  the  Inherent  generality  of  the 
results  for  TNT-llke  energy  fluxes.  Nevertheless,  for 
particular  accidents  and  containment  conditions,  the  rationale 
of  the  s low-release  results  can  frequently  be  employed  to 
yield  highly  defensible  containment  predictions. 
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9.10  Blast  Shielding  and  Inertial  Constraints.  All  of  the 
results  presented  In  Table  9.1  are  addressed  to  reactor  vessels 
with  no  external  constraints,  l.e.,  no  restraining  forces  other 
than  the  strain-generated  tensile  and  shearing  forces  In  the 
vessel  wall  Itself.  However,  In  the  design  of  many  reactors, 
there  exist  neutron  and  blast  shielding  of  various  types 
surrounding  the  reactor  vessel.  In  some  cases,  this  shielding 
and  other  external  constraints  can  enhance  excursion  containment 
potential  to  a  marked  degree.  In  the  EBR-II  reactor,  for 
example,  a  crushable  blast  shield  envelops  the  reactor  vessel; 
reference  (x)  Indexes  the  energy  absorption  potential  of  such 
a  shield.  Here,  radial  deformation  of  the  reactor  vessel  would 
be  opposed  by  the  crushable  shield,  and  the  opposing  force  would 
Increase  with  Increasing  vessel  strain.  This  type  of  external 
constraint  can  greatly  enhance  the  containment  potential  of 
reactor  vessels  for  the  entire  energy  flux  spectrum  of  nuclear 
accidents. 

Other  constraint  escaaples  are  the  graphite  neutron  shielding 
surrounding  the  Enrico  Fermi  reactor  and  the  large  pool  of  water 
surrounding  Oak  Ridge’s  HPIR  (High  Flux  Isotope  Reactor). 
Experiments  reported  In  reference  (l)  Indicate  that  the  graphite 
shielding  of  the  Fermi  reactor  enhances  Its  TOT  containment 
potential  by  a  factor  of  two  or  more.  The  same  enhancement 
would  result  from  the  HFIR  pool  water.  However,  these  masses  - 
the  graphite  and  pool  water  -  constitute  Inertial  constraints 
only;  hence,  they  are  effective  only  for  TNT  energy  fluxes  and 
those  of  like  order. 

No  general  statement  can  be  made  as  to  the  containment 
enhancement  factor  of  Inertial  constraints.  However,  Information 
presented  In  this  section  and  In  reference  (l)  can  be  employed  In 
some  cases  to  Increase  the  bounds  presented  In  Table  9.1. 
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FIG.  9.1  DEVIATIONS  OF  PREDICTED  CHARGE  WEIGHT  ( EQUATION  (9.1))  FROM  EXPERIMENTAL  RESULTS 


PENTOLITE  CHARGE  WEIGHT  (GM) 
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VESSEL  SIZE:  5"  ID  -  10"  LENGTH  -  1/8"  WALL 
VESSEL  MATERIAL:  304  STAINLESS  STEEL,  CAST 
CHARGE  LOCATION:  CENTROID  OF  VESSEL 
TEMPERATURE:  70°F 


FIG.  9.2  COMPARISON  OF  EXPLOSION  TESTS  IN  AIR -FILLED 
AND  WATER -FILLED  VESSELS 
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FIG.  9.3  VARIATION  OF  PRESSURE  AND  STRAIN  WITH  TIME  FOR  PROPELLANT-LOADED  VESSEL 


NOLTR  63-140 


VESSEL  SIZE:  5" ID -10"  LENGTH -1/4"  WALL 
VESSEL  MATERIAL:  304  STAINLESS  STEEL, CAST 


WATER  TEMPERATURE:  7(fF 


FIG.  9.4  COMPARISON  OF  PROPELLANT  AND  EXPLOSION  LOADINGS  IN 
WATER -FILLED  VESSELS 
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TABLE  9.1 

SAFE  *  CONTAINABLE  CHARGE  WEIGHTS 

FOR  REAL  VESSELS 

UNDER  VARIOUS  ACCIDENT  CONDITIONS 


QUANTITY  OF 
COOLANT 

IN  VESSEL 

n  ,  DURATION  OF  ENERGY  RELEASE  (milliseconds) 

(TNT),^  <  <  100  R.  /3 

D  in 

100  <  3n/R.  <  1000 

in 

SYSTEM  TEMPERATURE,  T  (°F) 

o 

rx 

H 

FULL 

NONE 

(Wj 

R  max 

1  <«R> 

T  =  200 

FULL 

NONE 

<®R>  ™x 

*2(W.) 

R  max 

*6  1  (W  )  1 

R  max  , 

*'•  1  <«R>  ™xr 

T  =  850 

FULL 

NONE 

(wj 

R  max 

*2  W„) 

R  max 

*3  1  ( W„  )  1 

1  R  '  maxi 

*3  1  (W  )  1+ 

1  R  max  1 

NOTES  1 .  Real  vessels  are  those  fitted  with  weldments,  nozzles, 
and  end  closures. 


2.  (TNT)[)  is  time  duration  af  TNT  energy  release. 

3.  R|n  is  the  numerical  value  of  R|(ft). 

4.  (WR)max  is  given  by  equation  (9.4). 

5.  l(W|^)^Qx  I  is  the  energy  density  (1050  calories  per  gram) 
af  (W|^)fy,ax  without  regard  ta  energy  flux. 

6.  *The  coefficients  af  these  quantities  are  hypothesized  values. 

7.  T  These  values  are  valid  only  far  the  case  of  extenjive 
metallic  heat  sinks  (see  section  9.8). 
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CHAPTER  10 

SUMMARY  AND  CONCLUSIONS  ^ 

10.1  Sximmary.  Nuclear  excursions  range  from  small  accidents 
to  thousands  of  megawatt-seconds  occtirrlng  In  times  from  a 
few  microseconds  to  seconds  or  more,  and  these  can  take  place 
anywhere  In  the  coolant-loss  domain  from  no  loss  to  a  complete 
loss.  The  violence  and  destructiveness  of  this  spectrum  of 
accidents  were  examined,  and  the  most  adverse  accident  (the 
most  difficult  to  contain  per  unit  energy  released)  was  found 
to  be  that  correlative  with  TNT  energy  flux  for  a  reactor 
vessel  completely  filled  with  liquid  coolant.  The  philosophy 
of  containment  and  Its  Increased  need  In  urban  areas  were  pre¬ 
sented:  from  these.  It  was  deduced  that  engineers  require  highly 
quantitative  Containment  Laws  to  effect  economical  reactor 
vessel  containment  design. 

The  principal  parameters  of  reactor  vessel  containment  were 
Identified  and  partially  characterized.  Consistent  with  Salnt- 
Venant's  semi-inverse  method  of  classical  mechanics,  a  synthesis 
of  the  lumped  parameters  was  employed  to  formulate  an  embryonic 
Containment  Law  format.  Extensive  experiments  In  Hopklnson 
models  (some  fabricated  with  weldments,  nozzles,  and  rigid  end 
closures)  of  reactor  vessels  were  conducted;  these  provided 
boundary  conditions  for  evaluating  the  embryonic  format  emd 
resulted  In  the  basic  TNT  Containment  Laws  presented  In  Chapter 
9.  The  consideration  of  less  adverse  accident  conditions  - 
energy  flxixes  less  than  that  of  TNT  and  loss  of  coolant  -  led 
to  hypotheses  of  Increased  and  more  reasonable  safe  upper  bo\mds 
for  accidents  of  this  type  (Table  9.1). 

10.2  Conclusions .  A  priori  -  fundamental  rationale,  extensive 
experimental  data,  and  comprehensive  syntheses  of  empirical  and 
theoretical  containment  knowledge  -  the  following  conclusions 
are  drawn: 
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(a)  The  basic  TNT  Containment  Laws,  equations  (9*1) 
through  (9.4),  and  their  attendant  conditions,  restrictions, 

and  qualifications  have  been  carefully  and  cautiously  generated. 
These  laws  are  known  to  be  consistent  with  extensive  experi¬ 
mental  results;  they  constitute  upper  boimds  that  are  wholly 
defensible  throughout  the  entire  spectrum  of  accidents  set 
forth  In  section  10.1. 

(b)  The  TNT  Containment  Laws  for  Real  Vessels  (those 
fitted  with  weldments,  nozzles,  and  rigid  end  closxircs), 
equations  (9.3)  and  (9.4), constitute  wholly  defensible  upper 
bounds  that  are  both  safe  and  reasonable  for  TNT-llke  energy 
flxixes  and  the  condition  of  no  liquid  coolant  loss  from  the 
reactor  vessel  prior  to  the  accident. 

(c)  The  mr  Containment  Laws  for  Real  Vessels,  equations 
(9.3)  and  (9.4), constitute  safe  upper  bounds  that  are  wholly 
defensible  throughout  the  entire  spectrum  of  accidents  set  forth 
In  Table  9.1.  The  reasonableness  of  those  bounds  decreases, 
however,  with  decreasing  energy  flux  and  Increasing  loss  of 
liquid  coolant  from  the  reactor  vessel  prior  to  the  accident. 

(d)  TNT  energy  flxix  and  no  liquid  coolant  loss  constitute 
the  most  adverse  accident  conditions.  Accidents  that  occur 
with  lesser  energy  fluxes  and/or  loss  of  liquid  coolant  ax*e 
less  severe  and  more  amenable  to  containment.  Increased  con¬ 
tainment  bounds  for  the  less  adverse  accident  conditions  are 
presented  In  Table  9.1.  These  bounds  are  hypothesized  values; 
they  are  believed  to  be  both  safe  and  x*easonable,  but  they  are 
not  as  highly  defensible  as  the  bounds  provided  by  the  TNT 
Containment  Laws  for  Real  Vessels,  equations  (9.3)  and  (9.4). 

(e)  The  Containment  Laws  and  other  bounds  presented  In 
this  report  az*e  expressed  In  terms  of  conventional  materials 
and  properties  -  quantities  readily  available  to  the  containment 
design  engineer. 
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(f)  Aside  from  reactors,  the  fundamental  rationale 

and  extensive  experimental  results  presented  herein  contribute 
directly  to  hlgh-straln-rate  phenomenology;  they  will  be 
applicable  to  pressure  vessels  generally  and  to  a  wide  spectrum 
of  dynamic  structural  response  problems. 

(g)  Since  the  TNT  (or  pentollte)  explosion  produces  shock 
and  blast  pressure  mechanisms  that  cover  a  wide  spectrvmi  of 
energy  flux,  and  since  the  fundamental  characteristics  of  these 
mechanisms  are  known,  their  applications  to  basic  excursion 
containment  studies  can  be  highly  profitable. 

(h)  The  reactor  containment  commvinlty  Is  seriously  enjoined 
and  cautioned:  In  general,  realistic  appraisals  of  the  excursion 
response  and  containment  potential  of  reactor  vessels  should  not 
be  attempted  without  basic  experimental  evidence  of  the  type 
presented  herein  and  In  reference  (l). 

The  basic  Containment  Laws  and  correlative  material  pre¬ 
sented  In  this  report  are  believed  to  constitute  the  most 
extensive  and  authoritative  reactor  vessel,  excursion  contain¬ 
ment  Information  existing  on  this  subject.  The  results  and 
conclusions  presented  In  this  report  are  final,  and  they  con¬ 
stitute  the  Laboratory's  official  opinion. 
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APPENDIX  A 

GENERAL  STRAIN  ENERGY  EQUATION 
FOR  EXCURSION -LOADED  REACTOR  VESSEL 


Consider  a  closed  right -cylindrical  vessel,  e.g.,  a 
reactor  pressure  vessel,  subjected  to  an  accidental  excur¬ 
sion  that  causes  the  vessel  wall  to  flow  rapidly  and  gross 
plastically  in  3-space  and  time  in  a  general  way.  If, 
during  the  dilation, 

(a)  the  pressure  vector  at  every  differential  surface 
of  the  internal  wall,  and 

(b)  the  strain  vector  of  every  differential  surface  of 
the  external  wall 

are  known,  what  is  the  strain  energy  correlative  with  contain¬ 
ment? 

We  shall  derive  the  general  strain  energy  equation  of 
dynamic  equilibrium  for  an  Internally- loaded,  hollow,  right- 
circular  cylinder  closed  with  rigid,  radial  constraints  at 
the  ends.  The  closures  are  right  discs  fitted  concentrically 
with  the  hollow  cylinder.  The  cylinder  and  closures  together 
are  defined  as  the  cylinder  unit  which  comprises  the  solid 
equilibrium  system  to  be  investigated.  The  cylinder  and 
closures  separately  constitute  the  primary  and  secondary 
systems,  respectively,  and  are  treated  as  phenomenological 
contlnua  for  which  it  is  assumed  only  that  the  material  is 
incompressible  and  of  uniform,  constant  mass  density;  the 
common  assumptions  of  homogeneity  and  Isotropy  will  not  be 
invoked  at  this  time.  The  dimensions  and  physical  properties 
of  the  secondary  system  are  such  that  it  is  not  stressed 
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beyond  the  elastic  limit  during  loading  of  the  cylinder  unit; 
thus  all  plastic  deformation  Is  restricted  to  the  primary 
system.  The  load  generating  and  transmission  media  (as  yet 
undefined)  internal  to  the  primary  and  secondary  systems 
comprise  the  fluid  or  tertiary  system. 


Throughout  the  entire  loading  It  is  postulated  that  all 
the  mechanical  surface-loading  work  done  on  the  solid  system 
manifests  Itself  in  the  Immediate  mechanical  forms  of  strain 
energy,  the  energy  required  to  accelerate  the  mass  of  the 
solid  system,  the  energy  required  to  compress  and  accelerate 
the  atmosphere  external  to  the  solid  system,  and  the  energy 
absorbed  or  expended  In  effecting  a  change  in  the  potential 
energy  of  position  of  the  solid  system.  For  any  time  t  let 


Epp  ...  be  total  surface  shearing-force  work  done  on 
solid  system  by  motion  of  pressure -loading 
forces  £arallel  to  Internal  wall,  ft-lb 


®pn  ...  be  total  surface  normal-force  work  done  on 
solid  system  by  motion  of  £ressure -loading 
forces  normal  to  Internal  wall,  ft-lb 

Ep  ....  be  total  £ressure-loadlng-force  work  done 
on  solid  system,  ft-lb 

El  ....  be  total  ^nertla-force  energy  generated  in 
solid  system  as  a  consequence  of  mass 
acceleration,  ft-lb 

EpQ  ...  be  total  change  In  potential  energy  of  position 
of  solid  system,  ft-lb 

Eg  ....  be  total  compressive  and  inertia-force  work 
done  on  £tmosphere  external  to  solid  system, 
ft-lb 


Eg  ....  be  total  £traln  energy  absorbed  by  solid  system 
as  a  consequence  of  deformation,  ft-lb 
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External  mechanical  work  done  upon  the  primary  and  secondary 
systems  must  be  Instantaneously  conserved  In  the  Immediate 
mechanical  forms  of  strain,  kinetic,  compression,  and  poten¬ 
tial  (position)  energy.  It  follows,  a  priori,  that  a  conse¬ 
quence  of  this  conservation  Is 

Ep  =  Epp  +  Epn  =  Es  +  Ei  +  Epo  +  Ea  (A-l) 

Plastic  straining  of  the  primary  system  Is  accompanied  by 
the  generation  of  heat  within  the  system  which  may  or  may 
not  be  dissipated  from  the  system  depending  upon  the  tempera¬ 
ture  of  Its  environment.  It  has  only  to  be  recognized  here 
that,  with  respect  to  equation  (A-l),  the  heat  of  elastic  and 
plastic  strain  Is  Inherently  accounted  for  In  the  strain 
energy  term  Eg.  Although  the  magnitudes  of  the  surface, 
strain.  Inertia,  compression,  and  position  energies  are 
functions  of  the  temperature  of  the  solid  system,  the  validity 
of  equation  (A-l)  is  Independent  of  temperature  and  heat 
flux  within  the  range  for  which  utilization  of  the  system  as 
a  pressure  vessel  remains  practicable . 

We  choose  r,  6,  and  z  to  be  the  radial,  circumferential, 
and  longitudinal  coordinates,  respectively,  for  a  cylindrical 
coordinate  system  In  which  z  Is  directed  vertically  upward. 

The  axes  of  revolution  of  the  cylinders  and  closures  coincide 
with  the  z-axls,  and  the  Internal  face  of  the  lower  closure 
lies  In  the  r,  0,  z  =  0  plane  as  shown  In  Figure  A-l.  Let  Tq, 
0Q,  and  Zq  be  Lagranglan  coordinates  which  denote  the  orig¬ 
inal  position  of  each  solid-system  particle  at  time 

t  =  t  Q  =  0  . 

A-l  Normal  Pres sure -Force  Work.  We  first  derive  an  expression 
for  the  normal  pressure-force  work  done  upon  the  solid 
system.  We  choose 
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r\  -  ri(0o,Zo  ,  t) 

0i  =  0i  (00  ,Zo  it)  (A-2) 

Zi  =  Zi(0o,  Zo  1  +) 

to  be  the  radial,  circumferential,  and  longitudinal  coordi¬ 
nates,  respectively,  which  define  the  internal  surface  of 
the  vessel .  These  functions  as  well  as  all  others  defined 
in  this  appendix  are  assumed  to  be  differentiable  to  the  extent 
required  for  the  work.  Let 


Pr"  Pr  ®  i 

,Hi  ,+) 

Pe=P®(n 

,  Hi  ,  +) 

(A-3) 

Pa=Pa(n 

,Hi  ,t) 

be  the  loading  pressures  generated,  respectively,  in  the 
r,  9,  z  directions  and 

d/\r=riddi  di\ 

d^e^dr\dz\  (A-4) 

dt^i  -  r,  dr\  d6\ 

be  the  projected  differential  areas  upon  which 

Pr  >  P®  1  Pa 

respectively,  bear  to  generate  the  differential  pressure 
forces 
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dP^  -  PrH  dQ\  dl\ 


o'Pe  =  Pe  dx\  di\ 

-  ^i'^\dr\  dQ\ 

Uj  =  Ujieo,  Zo,  t) 
Vi  =  Vil^o,  Ho,  t) 
Wj  =  Wj  (00,  Ho,  t) 


(A-5) 


(A-6) 


respectively,  define  the  radial,  circumferential,  and  longi¬ 
tudinal  displacements  of  the  solid  system  internal  surface, 
the  differential  displacements  through  which  the  pressure 
forces 


dP,,  dP,,  dP^ 


move  in  time  dt  are 


du 


d\i  ;--^  dQr,* 


dUi 


C^Ho+  • 


dUi 


i"  a0„  “‘'O’  az«“‘o’  at 


^  dv:  ^  dv:  ,  dVi  ^ 

dW:  dW:  dWj 

"'"1=  ar 


(A-7) 


We  elect  here  to  follow  the  path  of  a  single  particle  in 
time  so  that  the  Lagranglan  differentials 
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dQo, 


vanish,  and  we  obtain 


(A-8) 


The  differential  work  done  in  the  r,  9,  z  directions  is, 
respectively. 


dUj 


^^pnr'PrH  dd\  dZ\  d\ 


OV: 

c^Epne=Pa-^  dr\  di-^  d\ 


(A-9) 


dr\  dd\  d\ 


pn2  "  Pz  f'i 


It  will  become  evident  in  later  considerations  that  it  is 
highly  desirable  to  choose  differential  forms  which  will  per¬ 
mit  integrations  over  the  original  configuration  of  the  solid 
system.  To  this  end,  from  the  relations  (A-2),  we  express  the 
non-Lagrangian,  surface  coordinate  differentials 


dx\ ,  dB  \  ,  di  \ 


in  terms  of  Lagrangian  and  time  differentials  as 
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dd:  36:  36: 


^2i  -  .  ^2i 


3ii 


(A-IO) 


But  the  differential  areas 


dAf ,  dAg,  dP<i 


of  which 


dr\y  d6\,  di\ 


are  factorial  components,  are  functions  of  the  spatial 
coordinates 


r\ 


6 


I  > 


Zi 


alone,  so  that  the  time  differential  vanishes,  and  the  equa¬ 
tions  (A-IO)  reduce  to 


dr\  = 

d6\  = 

dl\  = 


3r-, 

3r\ 

CD 

o 

d6o+  ■ 

3io 

36 

-d6n  + 

36-^ 

36o 

W  W  Q 

dio 

3i\ 

■d6o  + 

dii 

dz, 


36 


die 


dir 


(A-11) 
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Combining  equations  (A-9)  and  (A-ll),  we  obtain 


^Epnr  =Pr  [  41^  44*"  C/Bq 


at 


a^o  <3^0 


/  a^i 

dz\ 

a^i 

dz\ 

azi  1 

(a.. 

^Ho 

+ - !- 

Ho 

a^o 

aH„ 

an! 

(A-12) 


^^pn©  “P©  O^t 


arj  azi 


3&q  3& t 


ddo  ddo 


dri  azi 


ddr 


awi 


/  ari 

\  a6„ 


aei 


an 

^Ho 

azi 

a^o 

(A-13) 

-o'! 

‘ari 

La^o 

ael 

ari 

azo 

—  o 

05  05 

ro  'O 

1 

an. 

(A-14) 

which  are  the  components  of  work  resulting  from  the  differen¬ 
tial  motion  of  an  internal  differential  area  of  the  solid 
system  subjected  to  dynamic  loading.  We  define 


O^Epn  “  O^Epnr 0^ Epne  O^E 


’  pnz 


(A-15) 


from  which 


■pn 


(^^pnr'*’ ^Epng  i/Epnj) 


't  •'A 


(A-16) 
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where  t  and’  k,  respectively.  Indicate  time  and  original 
Internal  area  limits  of  Integration.  It  Is  necessary  at 
this  point  to  note  that  the  change  In  variables  from  non- 
Lagranglan  to  Lagranglan  coordinates  represented  In  the 
Integral  (A-16)  Is  governed  by  the  general  transformation 


f(x,y)  dx.  dy  = 


a(x,y) 

<^(Xo,yo) 


dxo  dyo  , 


X  =  X  ( Xo^yo  ) 

y  =  y  (xo,yo) 


Writing  equation  (A-l6)  In  accord  with  the  above,  we  obtain 


auj/  dSi  aZj  de,  d^-,  \  .  aVi/aZi  dr-,  dz-,  ar, 


at\  de,  deJ^^^dAde^di; di,de. 


+  p/i 


awi  /  ar,  de-,  ar:  de-,  \~ 


* '  at 


dBo 


de^diod^ 


(A-17) 


which  Is  the  total  normal  pressure-loading-force  work  done 
on  the  solid  system. 

As  later  consideration  of  the  parallel  pressure-loadlng- 
force  work  Epp  requires  that 

EppE  0  (A-18) 

no  working  analytic  expression  will  be  written  for  It. 

A. 2  Inertia-Force  Energy.  We  next  derive  an  expression  for 
the  Inertia-force  energy  generated  In  the  solid  system.  We 
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define  p  to  be  the  uniform,  constant  mass  density  of  the 
solid  system,  and  the  mass  of  the  differential  volume  dV  Is 

c/m  =  p  dV  -  prQ  dx^  dB^  c/Zq 


If 


u  =u(ro  ,  @0.  Zo. 

v  =  v(ro  ,  @0  ,  Zo  .  (A-19) 

w  =  wlro  ,00  ,  Ho  1 


are  the  radial,  circumferential,  and  longitudinal  displace¬ 
ments,  respectively,  of  a  solid  system  particle,  the  corre¬ 
sponding  differential  displacements  through  which  the  mass  dm 
moves  In  time  dt  are 


C/U  = 


d\  - 


d^- 


0U 
d  Tq 

dv 

d  Xq 

d  w 

dro 


dro  + 

dxQ  + 

o'ro  + 


du 

d0o 

d  V 

d  w 


O'0o  + 

O'0o  + 

0^00  + 


du 

dzo 

dv 

(3ho 

d  w 

dZo 


O'Zo  + 

dio  + 

0'Zo+ 


du 

dt 

dv 

yr 

d  w 


O'! 

d^ 

d\ 


(A-20) 


We  elect  again  to  follow  the  path  of  a  single  particle  In 
time  so  that  the  Lagranglan  differentials 

dr^  ,  dB^  ,  diQ 


vanish,  and  the  particle  differential  displacements  and 
accelerations  may  be  written  as 
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du  = 

du 

at 

rfv-  rft, 

dv^  =  d1 

0\ 

(A-21) 

d^u 

d^u 

d\  _  a^  V 

d"w  _  a"w 

(A-22) 

■  ’ 

dt"  ^  at"  ’ 

c/t"  ^  at" 

Let 

Eir  »  Eie  »  Eiz 


be  the  translational  Inertia -force  energies  generated  by 
linear  acceleration  of  the  solid  system,  respectively.  In 
the  r,  0,  z  directions .  The  corresponding  differential 
quantities  may  be  expressed  In  accord  with  D'Alembert's 
principle  as 


rfEi,  =  rfro  rfSo  rfHo  rft 

p^Q  dxQ  cIBq  dio  d'\  (A-23) 

d^\i  =  yy  p^o  dxo  dd^  dzo  d1 

We  define 

fl'EiT=fl'Eir+  dE\g  +  C/E\t  (A-24) 


from  which 
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(fl^Ei|'  +  fl^EiQ+  ^Ejj) 


(A-25) 


Combining  eqiiatlons  (A-23)  and  (A-25)  and  simplifying,  yields 


•iT  - 


du  dw  \ 

IT / 


fl'Vo't 


(A-26) 


which  Is  the  total  translational  energy  of  the  mass  particles 
of  the  solid  system  or  the  total  energy  generated  as  a  conse¬ 
quence  of  linearly  accelerating  the  mass  particles . 

Let  X,  y,  and  Z  be  the  Instantaneous  principal  axes 
about  which  the  differential  mass  dm  rotates  and 


be  the  moments  of  Inertia  of  dm.  If 

a  =  a(ro,  00  ,  +  ) 

0  =  /3('b  ,  00  ,  Ho  ,  t )  (A-27) 

y  =  r(rQ  ,00  ,Zo,t) 

are  angular  displacements  about  X,  Y,  Z,  respectively,  the 
differential  angular  displacements  through  which  dm  moves  In 
time  dt  are 
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da  .  da  da  ,  da  , . 
da=  —  dro  +  —  c/Qq-^  -r- diQ-^  —  di 
dro  d9o  diQ  dt 


dy3  dy3 


ddQ  + 


dy-- 


dy 

d'(o 


i21 

d^o 


dXa^  -—dQa-^ 


d/3 

din  + 

d/3 

^Zo 

dt 

dy 

din  + 

dy 

dzo 

W  A.  Q 

dt 

d^ 


di 


(A-28) 


Consistent  with  the  Lagrangian  concept  of  motion,  the  spatial 
differentials  vanish,  and  the  angular  displacements  and  accel¬ 
erations  are  written  simply  as 


,  da  dP  ..  .  dy 

o'a  =  -^;^o'+,  o'/s  =-^ o'!,  o'y=-jp-o'+ 

d^a  _  d^a  d^B  _  d^$  d^y  _  d^y 
c/+"  "  d+"  ’  c/t"  ‘  dt"  ’  c/t"  ”  d+" 

Let 


(A-29) 

(A-30) 


^\X  >  ^\Y  >  EiZ 


be  the  rotational  Inertia-force  energies  generated  by  angular 
acceleration  of  the  solid  system  mass  particles  about  X,  Y, 
and  Z.  The  corresponding  differential  quantities  may  be 
expressed  In  accord  with  D'Alembert's  principle  as 


“  h 


d* a  da 


dt"  dt 


df 


dE\Y  -  ly 


d^B  dB 


dV  dt 


o'! 


(A-31) 
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dE\% 


at 


Determination  of  the  moments  of  inertia  appearing  in  equations 
(A"31)  requires  integration  over  the  differential  volxmie 
dV,  To  satisfy  this  requirement  it  will  suffice  to  treat  the 
”pie"  shaped  differential  volume  of  the  cylindrical  coordi¬ 
nate  system  as  a  rectangular  parallelepiped  with  sides 

Tq  d 9 q  ,  d^Q  ,  d  I.  Q 

The  moments  of  inertia  may  then  be  easily  expressed  as 


,  <yzc 
2 


2 


,  c/rp 
2 


dtp 

2 


^  tp  dQp 
2 


Tpd  Qp 
2 


(Y"+^")c^X  dX  d% 


Iv  --P 


U  --P 


,  O-Zo 

^  To  <3^00 

2 

r  ^  1 

f  ^ 

2 

2 

To 

2 

C^Ho 

^  dro 

2 

f  ‘ 

To  de^ 

2 

2 

2 

(X'+^Mo'X  d/Y  (A-32) 


(X"+  Y")  d/X  dY  d% 


Performing  the  Indicated  integrations,  the  equations  (A-32) 
yield 
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p  r 

Ix  =  ~\2  ''o  ^•'o  ddodlo  ( 

:o'ro)"+(o'Zo)"] 

Iy  =  if’  '^odro  0'5o^3'Ho| 

{x^dB^^  +  (o'Zo)"j 

p 

=  ,2  r^dx^dB^dt^ 

[^(o'ro)"  +  (roO'0o)"j 

We  define 


o'EiR=  o'EiY  +  dE\% 


from  which 


EiR  =1  I  (o'Eix  +  o'EiY  +  o'Ei^) 
't  •'V 


Oomblnlng  equations 
we  obtain 


EiR 


3^0  d£_ 

at"  at 

a"y  a  r 
at"  at 


(A-31),  (A-33),  and  (A-35)  and 


a"a  aa 

at"  at 


(o'ro)"  +  (o'Zo)" 


irodd^Y  +  {dloY 


]} 


{dro)^^{todeoY\  \dYd^  =  0 


(A-33) 


(A-34) 


(A-35) 

simplifying. 


(A-36) 
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which  is  the  total  rotational  energy  of  the  mass  particles 
of  the  solid  system  or  the  total  energy  generated  as  a  con¬ 
sequence  of  angularly  accelerating  the  mass  particles.  The 
total  energy  generated  in  the  solid  system  as  a  consequence 
of  mass  acceleration  is 


Ei  =  EiT  +  EiR=^ 


d^u  du  dv  dw  , 

at"  at  at"  at  at"  at  ' 


f 

I  . 

Ji  J' 


l'^  I  I  iu 

'V 


a"a  aa 
^t"  at 


(a'ro)2+(a'Ho)2 


a"/9  d0 


at  at  . 


(ro  ddoY^idirY 


+ 


a"y  ar 


at^  at 


[dxQ  Y  +  (ro 


} 


dV  o't 


(A-37) 


In  equations  (A-36)  and  (A-37)  each  component  of  rota¬ 
tional  energy  vanishes  as  a  consequence  of  the  first-order 
differentials  squared 

{dr^Y  ,  (ro  dB^^ ,  (o'Zo)" 
tending  to  zero  in  the  limit  as 

o'ro—  dB^-*  dZo—  0 

Important  physical  interpretations  can  be  imposed  upon  these 
squared  differentials .  It  is  postulated  that 

(a)  if  the  inertia  energy  generated  in  the  solid  system 
is  determined  by  analytic  integration,  as  In  equation  (A-26), 
all  of  the  energy  will  appear  in  the  form  of  translational 
energy,  the  rotational  energy  being  identically  zero. 
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(b)  if  the  Inertia  energy  generated  in  the  solid 
system  is  determined  by  nxxmerlcal  Integration,  the  squared 
differentials  do  not  tend  to  zero  but  must  be  assigned 
some  constant,  Incremental,  non-zero  value.  The  ratio  of 
rotational  Inertia-force  energy  to  translational  Inertia- 
force  energy  manifested  will  be  a  fvinctlon  of  the  dimensions 
assigned  to  the  elemental  voltme 


AV=  To  Aro  AOq  AZq 


As 


Aro—  AOq^  AZq—  0 


It  Is  clear  that 


Et 


0 


A-3  Potential  Energy  of  Position.  The  total  change  In  the 
potential  energy  of  position  of  the  solid  system  is 


E 


o'Vo't 


(A-38) 


A-4  Work  Done  on  Atmosphere ♦  To  complete  our  accounting  of 
the  energy  terms  appearing  In  equation  (A-l)  we  must  reflect 
upon  the  work  Eg  required  to  compress  and  accelerate  the 
atmosphere  external  to  the  solid  system.  We  first  recognize 
for  the  subject  experiments  that,  during  gross  motion  of  the 
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simulant,  pressure  acting  on  the  Internal  wall  will  always 
considerably  exceed  the  opposing  atmospheric  pressure.  We 
also  recognize  that  the  density  of  the  simulant  greatly 
exceeds  that  of  the  atmosphere.  Prom  elementary  gas  dynam¬ 
ics,  then,  we  can  deduce,  with  no  appreciable  loss  In  gener¬ 
ality  for  the  subject  work,  that  Eg  Is  a  negligibly  small 
quantity.  We  shall,  however,  retain  the  quantity  for  a 
later  and  more  appropriate  disposition. 

A-5  Exact  General  Equation.  Equation  (A-l)  may  be  written  as 

Es  ■  Ep  Ej  Epo~  Eg  =  Epp  +  Epn~E|j  Ejp  Epo~Ea  (A-39) 

which  states  that  the  strain  energy  absorbed  by  the  solid 
system  can  be  expressed  as  functions  of  the  dynamic  pressure 
bearing  upon  the  Internal  and  external  walls  of  the  system, 
the  mass  motion  of  the  system,  and  the  mass  motion  of  the 
atmosphere  external  to  the  system.  Combining  equations  (A-I7), 
(A-37),  (A-38),  and  (A-39)^  we  obtain 


Es  =  Epp  ■* 


duJ  dQj  dZ|  dQj  dZ| 


dl\de^  dz^  dz^  de^ 


de^dz^dl-P 


2 


(a-4o) 
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Equation  (A-4o)  Is  the  general  strain  energy  equation  of 
dynamic  equilibrium  for  a  vessel  subjected  to  Internal 
dynamic  loading. 

It  seems  In  order  at  this  time  to  make  a  further  state¬ 
ment  as  to  the  relation  between  the  general  strain  energy 
equation  and  correlative  experimental  work.  The  mathematical 
and  physical  formulation  of  the  strain  energy  equation  is  not 
a  theory  designed  for  experimental  confirmation;  but  rather 
the  principal  analytic  tool  through  which  experimental  work 
can  be  utilized  In  assessing  the  effectiveness  of  strain- 
energy  absorption  as  a  means  for  Increasing  the  containment 
potential  of  reactor  vessels.  This  assessment  can  be  expressed 
as  the  partition  of  mechanical  (strain  energy)  and  non-mechanical 
energy  resulting  from  excursion- type  energy  releases  In  simulated 
reactor  vessels.  With  the  total  energy  Input  known,  the  general 
strain  energy  equation  provides  the  means  for  obtaining  the 
desired  partition  without  further  knowledge  as  to  the  disposi¬ 
tion  of  non-mechanical  energy. 

A-6  Simplified  Working  Equations .  The  general  equation  (A-40) 

Is  greatly  simplified  for  the  case  of  the  homogeneous  and 
Isotropic,  right-circular  cylinder  closed  with  rigid,  radial 
constraints  at  the  ends.  The  conditions  of  circularity, 
cylindrlcallty,  and  rightness  of  the  vessel  provide  that  the 
spatial  limits  of  Integration  be  constant.  If  the  loading 
Is  rotatlonally  symmetric,  the  structural  response  of  the 
homogeneous,  isotropic,  solid  system  will  be  rotatlonally 
symmetric  with  the  result  that 


V  =  Pe=E 


PP 


and 


dQ\  di\  a/3  dy 


dzo  dOn  at  at 


=  0 


(A-4i) 
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dd] 

ddo 


(A-42) 


It  may  also  be  shown,  for  the  subject  work,  that  the 
change  in  the  energy  of  position  of  the  solid  system  is 
negligibly  small  relative  to  the  accompanying  generation  of 
pressure-force  and  inertia-force  energies.  With  respect  to 
the  total  energy  absorbed  at  any  time  t,  we  can  Impose  the 
condition 


P9 


c/ydi  =  0 


(A-43) 


with  a  loss  in  accuracy  of  less  than  one -thousandth  of  one 
per  cent.  We  also  elect  at  this  time  to  take 


Ea  =0 


(A-44) 


which  elementary  considerations  indicate  will  likewise 
result  in  a  loss  in  accuracy  not  greater  than  one-thousandth 
of  one  per  cent. 

Imposing  conditions  (A-4i),  (A-42),  (A-43),  and  (A-44) 
on  the  general  equation  (A-4o),  we  obtain 


dz,  dwj  dr\ 

-  +  D»i'i -  - 

dzo  df  dZo 


ddo  dZo  dt 


A-20 


NOLTR  63-140 


-P 


d  u  d^w 


To  dxo  ddo  dio  d-^ 


{dtoY  +  {di^^  To  dXo  d6o  dio  di 


(A-45) 


which  is  the  strain  energy  absorbed  by  a  right-circular, 
homogeneous,  isotropic  cylinder  subjected  to  internal,  rota- 
tionally  symmetric,  dynamic  loading.  To  express  explicitly 
the  integration  limits  of  equation  (A-45),  it  Is  necessairy  to 
write  distinct  integrations  over  the  primary  and  secondary 
systems.  We  are  particularly  Interested  in  the  case  for 
which  the  secondary  system  consists  of  right  discs  fitted 
concentrically  with  the  primary  system.  If  the  closures  are 
not  stressed  beyond  the  elastic  limit  during  loading  of  the 
cylinder  unit,  we  can,  with  a  negligible  loss  in  accuracy, 
assume  that  the  dimensions  of  the  discs  remain  constant  dur¬ 
ing  the  loading  phenomenon.  Imposing  upon  equation  (A-45)  the 
conditions  that  the  Internal  face  of  the  lower  closure  be 
constrained  to  lie  in  the  r,  9,  z  =  0  plane  and  that  the 
dimensions  of  both  closures  remain  constant,  we  obtain 


/•2Tr 


'  duj  dzj 


.^^]d9 

'dt  dZo/^®° 


diod^ 


’»i  /•Ho  /•r, 


•0 


-p 


d*u  du  d^w  dw  ,  ^  ^ 


'0  W  Ijo 
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2>»’/T,o 

/  d^a  da 

/  at"  at 

•'r,* 


(fl'ro)*+  (i/Zo)* 


ro  o^ro  c/eo  c/Zo  i/t 


+ 


P*''o  cfroddo  cff-PVi 


d‘w  dw 

dP’  at 


c/f 


(A-46) 


where 

IJ 

o  ....  Is  original  length  of  primary  system  between 
closures.  Inches 

rj^Q  ...  Is  original  Internal  radius  of  primary  system.  Inches 
1*60  •  •  •  original  external  radius  of  primary  system.  Inches 
to  ....  Is  time  at  t  »  to  =  0 

t]^  ....  Is  time  subsequent  to  t  =  to,  milliseconds 

V2  ....  Is  volume  of  upper  closure  of  secondary  system, 
cubic  Inches 

If  both  of  the  closures,  which  constitute  the  secondary 
system,  are  rigidly  constrained,  their  pressure  and  Inertia - 
force  energies  vanish,  we  take 


ill 


=  1 


(A-47) 


and  equation  (A-46)  reduces  to 


E,  = 


p*n 


awi  ari 


at  az( 


cf9„ 


C/Zo  Cff 
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3w  \ 

'dV' 


ro  dTo  dOo  dio  df 


da 

at 


{dtoV  +  {dio)*  I'd  ^^0  dlo 


(A-48) 


If  equation  (A-40)  is  integrated  over  the  entire  range 
"to  '^’a  "  (^1  )velocity * 0 


i.e.,  the  time  required  to  accelerate  and  decelerate  all 
particles  of  the  solid  system  from  and  to  a  state  of  rest, 
then  the  total  work  done  upon  the  system  in  accelerating  it 
must  be  precisely  equivalent  to  the  kinetic  energy  relin¬ 
quished  by  the  system  during  deceleration.  The  structural 
response  significance  of  this  statement  is  that  the  inertia 
forces  oppose  the  pressure  forces  and  decrease  the  Internal 
stresses  in  the  system  during  acceleration  and,  conversely, 
assist  the  pressure  forces  and  Increase  the  Internal  stresses 
during  deceleration.  The  net  effect  of  this  phenomenon  is 
such  that  with  respect  to  equation  (A-4o),  the  inertia  terms  Ej^ 
vanish.  It  follows  that  we  can  conveniently  reduce  equations 
(A-46)  and  (A-48),  respectively, to 


Esa" 


/  dui  dii  dwj  dri\ 

r'' at  aj„  '''''' at 


(A-4c, 
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for  the  case  of  constraining  the  lower  closure  and 


(3wi  dr\ 

PzH  ^1“  1 —  c/do  dZo  di 


df 


(3t  dir 


(A-50) 


for  the  case  of  constraining  both  closures. 

With  respect  to  the  general  energy  equation,  a  single 
important  simplifying  assumption  remains  to  be  considered. 

Let  the  loading  generator  be  a  circular  cylinder  with  length 
equal  to  and  longitudinal  axis  common  to  the  primary  system. 
Let  this  generator  be  capable  of  delivering  the  desired 
energy  flux  to  its  immediate  external  environment  In  a  rota- 
tlonally  symmetric  pattern.  If  the  annular  space  between  the 
external  wall  of  the  generator  and  the  Internal  walls  of  the 
primary  system  Is  filled  with  a  liquid,  it  is  assumed  for 
the  range 


millisecond  <  tg  i  100  milliseconds 

that  the  loading  is  delivered  normal  to  the  walls  of  the 
solid  system  and  that 


Pr  =  Ps=Pf  =  P  =  p(ri,  Zi  ,  +) 


(A-51) 


The  magnitude  of  the  error  incurred  by  this  assumption  has 
not  been  evaluated,  but  it  is  believed  to  be  small.  Imposing 
the  condition  (A-5l)  on  equations  (A-49)  and  (A-50),  resi)ectlvely, 
we  obtain 
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■so 


dzj  dwj 

dzo  ^  <5+ 


dQo  dZo  dt 


dfo  d9o  o't 


(A-52) 


for  the  case  of  completely  constraining  the  lower  closure, 
and 


■so 


d9o  dZo  dt 


(A-53) 


for  the  case  of  completely  constraining  both  closures.  Equa¬ 
tion  (A-52)  constitutes  the  basic  analytic  working  expression. 
It  may  be  seen  from  an  examination  of  this  equation  that  its 
solution  requires  experimental  determination  of 

(a)  the  pressure  bearing  upon  every  internal  differen¬ 
tial  area  of  the  solid  system  at  time  t,  l.e., 

p=p{ri  ,  2i  ,  t),  to<t<to  (A-54) 

(b)  the  radial  and  longitudinal  coordinates  of  every 
Internal  differential  area  of  the  solid  system  at 
t Ime  t ,  1 . e . , 


n  =  nl  2o  ,  +  ) 

Hi  =  Zi(  2oi 


,to<  t<  ta 


(A-55) 
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z 


FIG.  A  -  1  SECTIONAL  VIEW  OF  TYPICAL  REACTOR  SIMULANT 
DEFINING  COORDINATES  OF  ORIGINAL 
AND  DILATED  CONFIGURATIONS 
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APPENDIX  B 

STRAIN  PATTERNS  AROUND  AN  APERTURE  AND  NOZZLE 


This  appendix  presents  strain  data  that  provide  insight 
into  the  design  of  nozzle  Junctions.  The  deformed  vessels 
shown  in  Figure  7.6  were  accurately  measured  to  determine 
the  internal  radii  and  wall  thicknesses  in  the  regions  of 
discontinuity.  Figures  B-1  and  B-2  give  these  data  in  the 
form  of  graphical  profiles. 
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FIG.  B  -  1  WALL  PROFILES  OF  EXPLOSIVELY  DILATED  VESSEL  WITH  UNCONSTRAINED  APERTURE 
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FKJ.  8-2  WALL  P401US  Of  EXJ»LOSfVELY  DILATED  VESSEL  V^^fTH  SINGLE  NOZZLE 
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APPENDIX  C 

CALCULATION  OF  EMBRYONIC  FUNCTIONS,  AND  tg 


Table  C-1  presents  calculations  for  and  ♦2 
experiments  given  In  Table  7,1  for  which  containment  occurred. 
Nomenclature  and  units  are  compatible  with  equation  (8.I), 
and  tests  have  been  grouped  according  to  material.  Average 
mechanical  and  physical  properties  were  determined  for  each 
material  based  on  certified  standard  tensile  tests  (see 
section  7.1);  the  form  In  which  they  were  used  In  equation  (8,1) 
Is  presented  below. 


Type  Material 

Stress 

(^IJ 

®t 

Weight  Density 
(lb/ft3) 

304  stainless  steel 
(welded  plate) 

43,000  + 

151,700c 

490 

304  stainless  steel 
(centrlfugally  cast) 

30,000  + 

123,700c 

490 

4l0  stainless  steel 
annealed 

(centrlfugally  cast) 

66,000  + 

203,200c 

490 

4l0  stainless  steel 
hardened 

(centrlfugally  cast) 

151,500  + 

778,500c 

490 

212  flange  steel 
(welded  plate) 

46,500  + 

186,000c 

490 

1020  mild  steel 
(welded  plate) 

42,000  + 

150,000c 

490 

6061  -  T6  Alumlnvim 
(extruded) 

40,000  + 

78,000c 

170 

Table  C-2  presents  calculations  for  ?,  equation  (8.5). 


C-1 


CALCULATIONS  FOR  %  AND  FROM  EQUATION  ( 8. 1 ) 
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TABLE  C-1  (cont'd) 
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TABLE  C-2 

CALCULATIONS  FOR  ^  FROM  EQUATION  (8.5) 


TEST* 

NO. 

ID 

(in) 

WALL 

THICK. 

(  «n) 

CHARGE 

WEIGHT 

(gm) 

INITIAL 
STRAIN  RATE 
( in/in/ sec) 

e 

EQ.  (8.5) 

R  i/ho 
(ft/ft) 

DEVIATION  ** 
OF  STRAIN 

RATE  (%) 

304  STAINLESS  STEEL,  WELDED -PLATE  CYLINDERS 

2 

5 

0.060 

9.0 

2070 

3.34 

41.7 

+10.2 

9 

5 

0.118 

23.0 

2850 

2.27 

21.2 

+  0.4 

13 

5 

0.253 

56.1 

2900 

1.84 

9.9 

+  0.0 

18 

10 

0.118 

120.1 

1570 

3.14 

42.1i 

+  3.0 

22 

10 

0.253 

200.0 

1300 

2.45 

19. S 

+10.9 

25 

10 

0.253 

275.5 

1570 

2.51 

19.8 

+13.6 

26 

10 

0.253 

325.U 

1870 

2.36 

19.8 

+  6.8 

28 

10 

0.253 

359.8 

1940 

2.43 

19.8 

+10.0 

29 

10 

0.515 

669.3 

164o 

2.08 

9.7 

+13.7 

31 

10 

0.515 

775.5 

1970 

1.91 

9.7 

+  4.4 

33 

20 

0.253 

477.8 

568 

2.53 

39.5 

-14.0 

34 

20 

0.253 

960.0 

790 

2.90 

39.5 

-  1.4 

35 

20 

0.253 

960.8 

805 

2.85 

39.5 

-  3.1 

304  STAINLESS  STEEL,  CENTRIFUGALLY-CAST  CYLINDERS 


59 

10 

0.250 

189.8 

1450 

2.15 

20.0 

“  3.2 

60 

10 

0.250 

189.9 

1490 

2.09 

20.0 

-  5.9 

64 

10 

0.500 

481.2 

1810 

1.56 

10.0 

-15.2 

410  STAINLESS  STEEL,  CENTRIFUGALLY-CAST  CYLINDERS 


ANNEALED 

69 

5 

0.250 

28.0 

i860 

1.82 

10.0 

- 1.1' 

75 

10 

0.250 

160.6 

1310 

2.13 

20.0 

-  4.1 

81 

10 

0.500 

230.4 

950 

1.82 

10.0 

-  1.1 

212  FLANGE  STEEL,  WELDED  PLATE  CYLINDERS 


87 

10 

0.515 

229.3 

867 

1.92 

9.7 

+  4.9 

92 

20 

0.515 

1200.3 

670 

1.93 

19.4 

-11.9 

93 

20 

0.515 

1351.5 

650 

2.15 

19.4 

—  1.8 

MILD  STEEL,  WELDED  PLATE  CYLINDERS 

98 

LH 

1  0.253 

1  140.0  1 

1010  1 

1  2.48  1 

1  1 

1  +12.2 

NOTES 

*  Order  of  presentation  in  this  toble  is  some  as  Table  7.1. 

**  This  column  gives  deviotions  of  predicted  strain  rates  (equation  (8.7)  )  from  experimentol 
results  (column  5).  Positive  deviotions  indicate  predicted  strain  rates  greoter  than 
experimentol  volues;  conversely  for  negotive  deviotions. 
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APPENDIX  D 
HEAT-SINK  RATIONALE 

In  Appendix  0  of  reference  (l),  equations  were  developed 
for  a  hypothetical  model  that  prescribes  the  magnitude  of  the 
peak  Internal  blast  pressure  resulting  from  a  TNT  explosion. 
This  pressure,  as  defined  earlier  In  the  text.  Is  the  quasl- 
statlc  equilibrium  pressure  that  exists  In  a  closed  chamber 
after  decay  of  the  shock  wave.  The  reference  (l)  analysis  Is 
reproduced  here,  and  Its  relation  to  slower  release  rates  Is 
discussed. 

D-1  Hypothetical  Model.  Since  a  simple  quantitative  analysis 
of  an  explosion  In  a  closed  vessel  with  extensive  liquid  and/or 
metal  heat  sinks  Is  not  feasible,  the  following  h3rpothetlcal 
model  Is  entertained. 

(a)  A  closed  vessel  contains  a  charge,  heat  sinks,  and 
nitrogen  gas  with  the  following  Initial  conditions. 

Charge:  volume  Yq  and  temperature  T^  *  25°C 

Heat  sinks:  temperature  T^ 

Nitrogen:  pressure  Pq  *  14.7  psla,  volume  Vo»  and 

temperature  Tq 

(b)  The  charge  Is  detonated  In  a  manner  such  that  the 
product  gases  are  contained  within  the  solid  charge  volume  Vg, 
the  total  heat  of  detonation  Is  absorbed  and  stored  for  later 
distribution,  and  the  temperature  of  the  product  gases  becomes 
Tq  =  25°C.  The  total  heat  of  detonation  or  explosion  energy  Er 
momentarily  stored  Is 

Eip  -  Wh  (D-1) 


D-1 
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where 

W  ....  Is  charge  weight,  gm 
h  ....  Is  heat  of  detonation,  cal/gm 

Conditions  of  the  products  at  this  point  are  pressure  P^, 
volxime  Vc  temperature  T^,  and  they  are  related  by  the 
perfect  gas  law 


Pc  Vq  =  n^  RTq 


where 

n^  ....  Is  moles  of  gas  released  by  explosive 
R  ....  Is  universal  gas  constant 

(c)  Now,  part  of  the  explosion  energy  Ety  Is  given  back 
to  the  product  gases  to  raise  their  temperature  to  that  of  the 
surroundings  Tq,  or 

El  =  n^  Cy  (Tq  -  Tq)  (D-2) 


where 

El  ....  Is  energy  required  to  raise  product  gas  temperature 
to 

c^  ....  Is  specific  heat  at  constant  volvmie  of  charge 
products 

(d)  The  charge  products  and  the  Initial  nitrogen  gas 
are  mixed  In  an  Isothermal  process  for  which  the  mixed-gas 
pressure  Pi  Is 


P 


1 


(  ^  ) 


) 


D-2 


where 
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n2  ....  Is  stun  of  moles  of  Initial  N2  and  charge  product 
gas 

njj  ....  Is  moles  of  Initial  N2 

V2  ....  Is  sum  of  Initial  N2  volumes  and  charge  volvime 

(e)  A  portion  of  the  explosion  energy.  In  the  amount  E2, 
Is  now  released  as  a  shock  wave  which  deforms  the  reactor 
vessel  wall.  This  deformation  provides  an  additional  volxime 
V2  Into  which  the  mixed  gases  can  expand,  so  that  the  final 
volinne  Vp  becomes 


Vp  -  Vo  +  +  V2 


the  pressure  becomes 


Pi  (  ^  ) 


Vo 


(D-3) 


and  the  temperature  Is  still  Tq. 

(f)  Energy  Is  now  given  to  the  various  heat  sinks  In  the 
amount  E^.  Here  the  assinnptlon  Is  made  that  the  heat  sinks 
are  sufficiently  large  to  absorb  the  energy  E^  without  a  sig¬ 
nificant  rise  In  temperature.  Thus  the  mixed  gas  temperature 
still  remains  essentially  at  Tq. 

(g)  The  unused  explosion  energy  Ep  at  this  point  Is 

Ep  =  Eqi  -  E2  -  E2  - 


D-3 


(D-4) 
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This  remaining  energy  is  given  to  the  mixed  gases  to  raise 
the  gas  pressure  to  a  new  value  Pp  (peak  Internal  blast 
pressure).  If  the  compressibility  of  the  heat  sinks  Is 
neglected,  then 


Er  (k  -  1) 

Vp 


+  P2 


where  k  Is  the  ratio  of  specific  heat  at  constant  pressure 
to  the  specific  heat  at  constant  volvmie.  Substitution  of 
equations  (D-l),  (D-2),  (D-3),  and  (D-4)  Into  this  expression 
yields 


ni 


ncR 


"o  (  4  )  (  ^  )  -  ^  (To  -  T,) 


+  r  (Wh  -  Eg  -  E3) 


(D-5) 


Equation  (D-5)  prescribes  the  peak  Internal  blast  pressure  for 
the  hypothetical  model. 

D-2  Special  Cases.  If  It  is  assumed  that  the  liquid  and/or 
metal  heat  sinks  are  sufficiently  large  to  absorb  all  of  the 
energy  released  by  the  explosive,  then 

(Wh  -  Eo  -  Eo)  -  —  (Tq  -  T.)  =  0 
Vp  '  3^  Vp  '  o 
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and  (D-5)  reduces  to 


P 


( 


m 

"N 


)  (  —  ) 

’  '  Vp  ^ 


(D-6) 


Experimental  results  with  the  water,  molten  sodium,  and/or 
metal  heat  sink  models  presented  In  reference  (l)  exhibit  good 
agreement  with  equation  (D-6). 

For  marginal  containment  In  reactor  vessels,  the  charge 
volinne  Vq  Is  usually  very  small  compared  with  the  sum  of  the 
Initial  gas  volume  Vq  and  the  deformation  volume  V^;  for  this 
case 


Vp  «  Vq  +  Vg 

Conversely,  the  gas  released  by  the  charge  n^.  Is  usually  very 
large  compared  with  the  amount  of  Initial  gas  njj;  this  yields 


ni  «  n^ 

In  view  of  the  above  simplifications,  equation  (D-6)  can  be 
approximated  by 


P 


P 


>  (  ^  ) 


Vo 
Vo  + 


(D-7) 


Since  n^j  Is  proportional  to  the  explosive  charge  weight  W, 
and  since  the  perfect  gas  law  gives 


-  RTo 

Tiff  o 
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equation  (D-7)>  for  a  given  vessel,  can  be  expressed  as 


Pp  «  W  (D-8) 

D-3  Slower  Release  Rates.  The  above  analysis  was  based  on 
a  TMT  explosion.  Entertain  now  a  release  rate  slower  than 
that  of  TNT.  The  assumptions  are  that  the  non-explosive 
charge  weight  W  would  be  sufficient  to  produce  marginal  con¬ 
tainment,  that  the  heat  sinks  would  be  sufficiently  large  to 
absorb  energy  without  a  significant  rise  In  temperature,  emd 
that 


^c  <  <  Vo  +  ^2 


Although  most  reactors  possess  the  extensive  metal  and  coolant 
heat  sinks  In  the  region  of  the  core  necessary  to  satisfy  the 
above  assiamptlon  on  energy  absorption.  It  Is  shown  In  reference 
(l)  that  liquid  coolant  or  metal  heat  sinks  separately  may  not 
be  sufficient. 

For  marginal  containment,  the  vessel  deformation  pattern 
of  slow  releases  would  be  similar  to  that  produced  by  TNT  In 
a  liquid  filled  vessel;  the  deformation  volumes  V2  would  be 
equivalent.  The  slower  release  would  have  to  be  accompanied 
by  energy  In  the  amount  E2  to  produce  this  deformation.  The 
event  time  would  be  longer  for  the  slower  rate  and  would  allow 
more  time  for  heat  transfer.  Since  In  the  short  event  time  of 
the  TNT  release  all  of  the  remaining  energy  Is  transferred  to 
the  heat  sinks,  certainly  no  less  could  be  expected  for  the 
slower  release. 


D-6 
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Prom  an  examination  of  the  above  conditions  in  light  of 
the  preceding  analysis,  it  is  hypothesized  that  equation  (D-8) 
adequately  prescribes  the  equilibrium  pressure  in  a  reactor 
vessel  with  extensive  heat  sinks  irrespective  of  the  energy 
release  rate. 
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